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We study the formation and destruction of molecules in the ejecta of Popu- 
CN ! lation III supernovae (SNe) using a chemical kinetic approach to follow the evo- 

, , lution of molecular abundances from day 100 to day 1000 after explosion. The 

chemical species included in the study range from simple di-atomic molecules to 
more complex dust precursor species. All relevant molecule formation and de- 
struction processes that are unique to the SN environment are considered. Our 
q work focuses on zero-metallicity progenitors with masses of 20, 170, and 270 M , 

£3 j and we study the effect of different levels of heavy element mixing and the in- 

c3 ', ward diffusion of hydrogen and helium on the ejecta chemistry. We show that 

the ejecta chemistry does not reach a steady state within the relevant timespan 
(~ 3 yr) for molecule formation, thus invalidating previous results relying on 
t— 1 I this assumption. The primary species formed in the harsh SN environment are 

^ I O2, CO, SiS, and SO. The SiO, formed as early as 200 days after explosion, is 

rapidly depleted by the formation of silica molecular precursors in the ejecta. 
The rapid conversion of CO to C2 and its thermal fractionation at temperatures 
G\ I above 5000 K allow for the formation of carbon chains in the oxygen-rich zone of 

the unmixed models, providing an important pathway for the formation of car- 
I bon dust in hot environments where the C/O ratio is less than 1. We show that 

the fully-mixed ejecta of a 170 M progenitor synthesizes 11.3 M of molecules 
whereas 20 M and 270 M progenitors produce 0.78, and 3.2 M of molecules, 
respectively. The admixing of 10% of hydrogen into the fully-mixed ejecta of 
the 170 Mq progenitor increases its molecular yield to ~ 47 M . The unmixed 
ejecta of a 170 M progenitor supernova without hydrogen penetration synthe- 
sizes ~ 37 M & of molecules, whereas its 20 M counterpart produces ~ 1.2 M & . 
This smaller efficiency at forming molecules is due to the large fraction of He + in 
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the outer mass zone of the ejecta. Finally, we discuss the cosmological implication 
of molecule formation by Pop. Ill SNe in the early universe. 

Subject headings: astrochemistry — supernovae: general — early universe - 
molecular processes 

1. INTRODUCTION 

Large column densities of dust are required to explain the reddening of background 
quasars and Lyman a systems at high redshift (z > 6) (Pettini et al. 1994, Pei & Fall 1995) 
and about 2 x 10 8 M Q of dust is derived from the infrared (IR) spectrum of the hyperluminous 
galaxy SDSS Jl 148+5251 at redshift z = 6.4 (Bertoldi et al. 2003, Robson et al. 2004, Beelen 
et al. 2006, Dwek et al. 2007). The origin of such large quantities of dust when the universe 
was less than 1 Gyr-old is still a matter of debate. In the local universe, dust forms in 
high density and temperature regions encountered in circumstellar environments such as the 
winds of Asymtotic Giant Branch (AGB) stars and supergiant stars, the colliding winds 
of Wolf-Rayet stars, and finally, the ejecta of core-collapse supernovae (CCSNe). In our 
Galaxy, most of the dust is produced by low-mass stars ascending the AGB. However, their 
long evolutionary main-sequence lifetime (a few Gyrs) excludes them from being possible 
dust contributors at high redshift. Conversely, very massive stars evolve much more rapidly 
(time scales ~ 1 Myr), and can be possible dust makers in the early universe. As to the first 
generation of stars, hereafter Pop III stars, they are expected to be very massive (Omukai 
& Nishi 1998, Abel et al 2002, Bromm et al. 2002). Indeed in the absence of metals, the 
cooling in primordial clouds is only provided by molecular hydrogen and thus precludes 
efficient gas fractionation. These Pop. Ill stars firstly need to synthesize heavy elements by 
thermonuclear reactions in their cores and reach their explosive ends to possibly condense 
dust in their massive supernova ejecta. Therefore, pair-instability supernovae (PISNe) are 
perhaps the first dust contributors to the pristine, young universe. 

The build-up of a molecular phase is a pre-requisite to dust nucleation and condensation. 
Indeed, it provides the molecular precursors from which dust forms and its composition 
depends on the initial elemental composition of the gas and the physical processes pertaining 
to it. Furthermore, molecules produced in the early universe can have an important effect on 
the cooling of the interstellar medium. Molecules have been detected in low-redshift CCSNe 
as early as 100 days post-explosion. Specifically, the IR ro-vibrationals transitions of CO 
and SiO were detected in SN1987A (Catchpole & Glass 1987, Meikle et al. 1989, Roche et 
al. 1991), CO fundamental bands were observed in the Type II SNe SN1995ad (Spyromilio 
& Leibundgut 1996), SN1998s (Gerardy et al. 2000) and SN202dh (Pozzo et al. 2006) whilst 
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SiO detection was reported by Kotak et al. (2006) in SN2005af. More recently, detection of 
CO with the Spitzer satellite in Cas A, a 300 year-old supernova remnant, is reported by 
Rho et al. (2009). It is therefore reasonable to expect molecules to form in the ejecta of 
massive Pop. Ill star supernovae. 

Existing models for CO and SiO formation in SN1987A ejecta only considered a limited 
number of chemical processes applicable to low-temperature gases and assumed that steady 
state held for chemistry in calculating the evolution of molecular masses (Petuchowski et 
al. 1989, Lepp et al. 1990, Liu & Dalgarno, 1994, 1995, 1996, Clayton et al. 1999, 2001, 
Geahard et al. 1999). Models for the formation of dust in SN1987A and PISNe used a 
classical nucleation theory approach to describe the growth of solids in the ejecta, ignoring 
the nucleation stage of dust, in which gas phase species and dust molecular precursors are 
formed (Kozasa et al. 1989, Toddini & Ferrara 2001, Nozawa et al. 2003 (hereafter NK03), 
Schneider et al. 2004 (hereafter SFS04)). Since the amount of carbon and oxygen locked 
up in CO is a major factor in determining the dust composition, Toddini & Ferrrara (2001) 
and SFS04 do consider the formation of CO in fully microscopically-mixed ejecta. However, 
their treatment is extremely over-simplified as their CO chemistry includes only two chemical 
processes and is assumed once again at steady state. In a first attempt to model molecular 
formation with a chemical kinetic approach in fully microscopically-mixed Pop. Ill SN ejecta, 
Cherchneff Sz Lilly (2008) (herafter CL08) show that the chemistry does not reach a steady 
state over the timespan studied. Furthermore, they identified other chemical processes than 
those considered by SFS04 that are of paramount importance to the formation of CO and 
other molecular species in the ejecta. 

In this paper, we study the formation and evolution of molecules, including gas-phase 
dust precursors, in the ejecta of Pop III SNe. We define a large set of chemical reactions 
relevant to the dense and hot SN environment but applicable to other circumstellar media 
as well. In addition to this extensive reaction network, with new updated rates compared 
to those used by CL08, we include processes unique to the radioactive environment of SN 
ejecta: destructive processed induced by Compton electrons created by the down-scattering 
of 7-rays, and by ultraviolet (UV) photons emitted by collisionally-excited atoms and ions. 
We consider fully microscopically-mixed and unmixed ejecta for different progenitor masses. 
In addition, the presence of hydrogen and helium can dramatically affect the chemistry of 
the ejecta. We therefore examine the effect of the inward diffusion of hydrogen and its effect 
on the molecular yield of the SNe. 

The molecular budgets of SN ejecta are determined by their physical conditions and 
composition. In §2 we describe the relevant parameters: explosion energy, initial density, 
the ejecta mass, and composition of the supernovae under study. These are used to fol- 
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low the evolution of ejecta temperature and density as a function of time. The reactions 
taking place in the different layers of the ejecta depend on the degree of elemental mixing, 
and the section also describes the prescriptions we used to characterize this process. The 
different mechanisms pertaining to the ejecta chemistry are described in §3. We first give 
a brief overview of the mathematical formalism, followed by a detailed description of the 
non-thermal destructive processes that operate in a SN environment. The results of our cal- 
culations for mixed and unmixed ejecta of various mass progenitors are presented in §4, and 
in §5 we briefly summarize and discuss the different implications of these results. Models for 
the formation of dust precursors in similar environments will be presented in a forthcoming 
paper (Cherchneff & Dwek in preparation). 



2. MODEL INPUT PARAMETERS 

2.1. Physical Conditions of the Ejecta 

Currently, there are no observational constraints on the Pop III supernovae events and 
the evolution of their explosive ejcta. We therefore base our ejecta models on available 
theoretical explosion models, in particular those of Heger & Woosley (2002), Umeda & 
Nomoto (2002, hereafter, UN02) and NK03. Using these models, we derive simple analytical 
expressions for the gas parameters such as temperature, number density and velocity and 
consider various levels of mixing in the ejecta. Three progenitor masses are studied: two 
massive progenitors of mass 170 M Q and 270 M Q chosen as surrogates to PISNe, and one 
low mass 20 M Q progenitor chosen to describe primordial CCSNe. 

The gas temperature T in the post-explosion gas is determined mainly by the explosion 
energy which is released as kinetic energy into the gas. NK03 present various models of ejecta 
for CCSNe, PISNe and hypernovae. Each model is characterized by an explosion energy, a 
mass cut marking the division between the matter which remains in the core and that which 
is ejected and the 56 Ni mass produced in the explosion above the mass cut. By solving 
the radiative transfer equation taking into account the energy deposition by radioactive 
elements, NK03 derive temperature profiles for each of their massive progenitor ejecta. For 
our PISNe models, we chose the NK03 temperature profile of the oxygen-rich region for their 
170 M Q PISN unmixed case to describe our 170 M Q ejecta temperature variation with time. 
The temperature characterizing the helium core of PISNe modeled by Fryer et al. (2001) 
stays almost constant over the core extent. We then assume a constant temperature in the 
inner He core region so that temperature is independent of mass coordinate M r and simply 
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fit the profile by the following power law of degree five 

T(t) = T x (A- Bx + Cx 2 - Dx 3 + Ex 4 - Fx 5 ), (1) 

where To = 21,000 K is the gas temperature at some fiducial time to = 100 days, x = t/to, 
and the fitting coefficients A, B, C, D, E and F are equal to 1.699, 8.568 x 10 -1 , 1.761 x 10 -1 , 
1.762 x 10~ 2 , 8.229 x 10" 4 , and 1.331 x 10~ 5 , respectively. Our 270 M temperature profile 
is assumed to follow the same variation with time as given in eq. (pQ). To account for the 
greater kinetic energy imparted by the explosion of a 270 M progenitor, we multiply To by 
a multiplying factor of 1.5. This value corresponds to the ratio of the central temperatures 
for the 170 M and 270 M & progenitor models of Heger & Woosley (2002). 

For our CCSNe model, we assume that the gas temperature follows the variation of the 
C20 unmixed model of NK03. Assuming the ejecta follows quasi-adiabatic expansion, the 
temperature evolution with time is given by 

/ t \ 3(1-7) 

T(t) = To x (-) , (2) 

where 7 is an 'adiabatic' index. Using eq. (J2J) to fit the NK03 C20 temperature profile gives 
a 7 value of 1.593. 

The ejecta expansion becomes homologous a few hours after explosion so that the gas 
density varies with time according to 

p(M r , t) = p{M r , t ) x (£j , (3) 

where M r is the mass coordinate. As for temperature, and according to the PISN models 
of Fryer et al. (2001), we assume a constant gas density over the helium core so that no 
dependence with the mass coordinate is considered. The gas density profile chosen in the 
present study is that of NK03 for their 20 M and 170 M progenitors, assuming p(600 days) 
equals to 3 x 10 -14 g cm -3 for all progenitor masses we study. 

The gas number density n(t) for each model is given by: 

n(Mr, t) = p(M r , t)/p gas (t) (4) 

where p gas {t) is the mean molecular weight of the gas, which varies with time as the chemical 
composition of the ejecta changes due to molecule formation. 

For the sake of simplicity, we define a constant ejecta velocity v determined by the 
explosion energy of the progenitor through its conversion to kinetic energy and given by 
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where Eq is the explosion energy and M e j is the mass ejected above mass cut during explosion. 
Explosion energies for the present models are listed in Table 1. The respective ejected masses 
are assumed to be equal to the progenitor masses as Pop. Ill stars do not experience mass loss 
during their evolution due to the lack of dredged-up metals in their photospheric composition 
and the consequent wind acceleration through metallic lines. The velocity is assumed to be 
constant over the mass zones in the ejecta. 

We are interested in studying the chemistry from time to = 100 days to t = 1000 days. 
The initial time range is justified by the appearance of CO, SiO and dust as early as 110, 160, 
and 450 days, respectively, in the ejecta of SN1987A (Catchpole & Glass 1987, Danziger et 
al. 1991, Wooden et al. 1993). Our PISN ejecta are hotter than that of SN1987A but similar 
temperature regimes are encountered at times t > 300 days. The final time is determined by 
the time when the gas density and temperature in the ejecta are too low to foster efficient 
gas-phase molecular formation. The ejecta parameters are summarized in Table 1 whereas 
the variation of the gas density and temperature with time t is illustrated in Figure 1 for 
the various progenitor masses considered. 



2.2. Mixing in the Ejecta 

Mixing is likely to occur during explosion for the instability of the nickel bubble resulting 
in the development of Rayleigh- Taylor instabilities (Woosley 1988, Arnett 1988, Herant & 
Benz 1991, Miiller et al. 1991, Kifonidis et al. 2003) and instabilities developing in the 
post-shocked regions of the propagating blast wave (Chevalier 1976, Bandiera 1984, Benz & 
Thielemann 1990). Evidence for strong 56 Ni mixing in the shell of SN1987A is found through 
observation of hard X-rays stemming from 56 Ni/ 56 Co decay and 7-rays down-scattering as 
early as 140 days after explosion (Itoh et al. 1987, Pinto & Woosley, 1988, Sunyaev et 
al. 1990). Another direct evidence for heavy element mixing comes from the extraction 
and study of isotopically anomalous inclusions in meteorites. The presence of type X silicon 
carbide (SiC) grains and silicon nitride (S3N4) inclusions bearing the 44 Ti supernova signature 
in meteorites suggests deep and inhomogeneous mixing between the various heavy element 
mass shells (Zinner 2006). As to light elements, hydrogen deep mixing down to the inner 
mass zones is invoked to explain the plateau shape of SN1987A light curve at times greater 
than 80 days (Woosley 1988, Arnett & Fu 1989, Shigeyama & Nomoto 1990). However, H 
mixing is likely to occur at a macroscopic level with the formation of H-rich bubbles in the 
unstable layer located between the He-CO and the H-He interfaces as early as a few hours 
after explosion (Fryxell et al. 1991, Herant & Benz 1991, Miiller et al. 1991, Burrows & van 
Riper 1995, Kifonidis et al. 2003). Some hydrogen diffusion may also occur at the base of the 
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hydrogen envelope where helium and heavy element-rich fingers are simultaneously present 
(Miiller et al. 1991, Kifonidis et al. 2003). To circumvent the problem of the complexity 
of mixing, NK03 considered explosive models with and without mixing, which result in two 
extreme cases for their ejecta: a fully mixed gas and a stratified ejecta in which each layer 
reflects the prior nucleosynthesis stages of the progenitor, except for the inner most layer 
which is the result of explosive nucleosynthesis. We follow the same strategy in the present 
paper: for the two very massive progenitors, 170 M and 270 M and the core-collapse SN 
progenitor, 20 M , we consider fully mixed ejecta while stratified ejecta are studied for the 
170 M and 20 M Q progenitors. 



2.3. Initial Ejecta Composition 

The initial chemical compositions for the primordial supernova ejecta models are taken 
from fully mixed and unmixed explosion models. For fully mixed ejecta, mass yields are from 
UN02 and Heger & Woosley (2002) for our massive 170 M & and 270 M & progenitors and 
from UN02 for the 20 M & CCSN model. The unmixed ejecta compositions for the 170 M Q 
progenitor and the 20 M progenitor are those of NK03. We subdivide their helium cores in 
zones of distinct chemical composition. For the 170 M progenitor, we consider five distinct 
zones: (1) is Si/S/Fe-rich from to 20 M Q , (2) is O/Si/S-rich from 20 to 40 M , (3) is 
O/Mg/Si-rich from 40 to 55 M , (4) is O/C/Mg-rich from 55 to 78 M , and finally (5) is 
O/C/He-rich from 78 to 82 M & . For the 20 M progenitor, the zoning comprises four zones 
as follows: (1) is Si/S/Fe-rich and extends from 2.4 to 3 M , (2) is O/Si/S-rich from 3 to 
3.6 M , (3) is O/C/Mg-rich from 3.6 to 4.95 M , and finally (4) is He/C/O-rich from 4.95 
to 5.85 Mq. For the sake of simplicity, we ignore the variations of the elemental abundances 
with mass coordinates and take them to be constant within each zone. 

We use these elemental mass yields to calculate the total number of elemental species, 
the gas initial molecular weight {i gas (t ) and gas number density n(M r ,t ). We then derive 
the number density at 100 days post-explosion for each element, and use these data as initial 
conditions when solving our set of non-linear, coupled, differential equations (see section 3). 
These data along with the resulting initial gas mean molecular weight are tabulated in Table 
4 for fully-mixed ejecta whilst Table 5 summarizes the post-explosion chemical composition 
as a function of zoning for our unmixed ejecta models. We then explore the effect of hydrogen 
mixing by diffusion from the progenitor envelope for the 170 M Q fully-mixed ejecta. In doing 
so, we set the H content as a free parameter and assume that H can microscopically mix 
within the heavy element-rich He core. Values of H mixing ranges from 0% to 10% of the 
total hydrogen envelope mass given by UNO 2. 
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3. 



EJECTA CHEMISTRY 



A chemical kinetic description of the ejecta is based on a gas initial chemical composition 
and a set of chemical reactions describing the chemical processes at play and applied to the 
ejecta physical conditions. For the high gas temperatures and densities characterizing our 
modeled ejecta, they include: (1) termolecular reactions efficient in high density media and 
where formation of molecules occurs through collision with the ambient gas which carries 
away the excess energy of the reaction. We also consider their reverse processes which are 
thermal fragmentation through collisions with the ambient gas at high temperatures; (2) bi- 
molecular processes, predominantly neutral-neutral (hereafter NN) reactions. The reactions 
with activation energy barriers require high temperatures to overcome the energy barrier 
whereas reactions without activation energy can proceed at lower temperatures; (3) ion- 
molecule reactions (formation/destruction reactions and charge exchange processes), which 
have no energy barrier and can therefore contribute at low gas temperatures; (4) temperature- 
independent radiative association reactions (hereafter RA) in which the formation of a species 
occurs through the emission of a photon which carries off the excess energy released in the 
formation of the adduct. 

Because of the unique nature of SN ejecta, being powered by the decay of radioactive 
elements and its high temperatures, we also consider the effects of various non-thermal 
processes including the destruction of molecules by the cascade of energetic electrons and 
UV photons that are generated by the down-scattering of 7 rays in the ejecta. These non- 
thermal processes will be discussed in more details in the following section. 

The temporal variation of the number density of a molecular specie % located in a given 
mass zone M r is described by the following rate equation: 



where P; are the production (= formation) processes and L x are the loss (= destruction) 
processes for species i, and k^ and k^ are the temperature-dependent rates for reactions 
between species i-j and i-k, respectively. The reaction rates kji and k ik are written under 
the form of Arrhenius expressions such as 



where T is the temperature given by equations ([I]) and ([2]), Ay the Arrhenius coefficient in 
s -1 molecule -1 , cm 3 or cm 6 s _1 molecule -1 for a unimolecular, bimolecular or termolecular 
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processes respectively, v reflects the temperature dependance of the reaction, and is the 
activation energy barrier in K _1 . The ensemble of equations as in eq.([6]) represents a set of 
N non-linear, coupled, ordinary differential equations to solve, where N is the number of 
species (atoms, ions, molecules, and electrons) included in the chemical description of the 
ejecta. In total, the system comprises up to 79 species listed in Table 2 and between 400 
to 500 reactions, depending on the ejecta region under study. The full chemical scheme is 
listed in Table 9. Reaction rates have either been measured under laboratory conditions or 
theoretically calculated using transition state theory. When not documented, the rates are 
estimated using 'educated' guesses. The NIST chemical kinetics database is used as primary 
source for NN processes, completed by the UDFA06 database (Woodall et al. 2007) when 
necessary. The numbering of some reactions specified in the text refers to that of Table 9. 
For dust precursor formation, processes and reaction rates stem from combustion chemistry 
studies, environmental and material sciences, and full details are given in Cherchneff & Dwek 
(in preparation). 



3.1. Non-thermal Processes 

3.1.1. Destruction by Compton Electrons 

A supernova explosion produces 56 Ni which rapidly decays to 56 Co with a half life of 
6 days, which in turn decays to 56 Fe with a half-life of ~ 77 days. The decay of 56 Co 
deposits 3.57 MeV in the form of 7-rays in the ejecta, which powers the SN light curve. 
(Woosley et al. 1989). Compton scattering degrades the 7-rays, which have an average 
energy of 1.24 MeV, to hard X-rays which through a cascade of inverse Compton, ionization 
and recombination processes degrade further into UV photons. The fast Compton electrons 
thermalize by heating, exciting and ionizing the ejecta, adding to the reservoir of UV photons. 
In SN1987A , the light curve between days 100 and 1000 could be reproduced if 0.075 M G of 
56 Co was ejected in the explosion. 

The fast Compton electrons and UV radiation can have a significant effect on the chem- 
istry of the ejecta. Compton electrons were proposed in several studies to be one of the 
dominant destruction routes to molecules in SN ejecta (Liu & Dalgarno 1994, 1995, 1996). 
Clayton et al. (1999, 2001) furthermore suggested that atomic carbon and carbon dust can 
be produced in oxygen-rich part of the ejecta by CO dissociation due to collision with Comp- 
ton electrons. We therefore pay special attention to deriving destruction rates for similar 
radioactivity-induced processes in the ejecta of our primordial SNe. In addition, we will also 
explore the role of UV photons on the chemistry of the ejecta. 
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Not all the radioactive decay is deposited in the ejecta. The early emergence of 7- 
rays and X-rays a few hundred days after the explosion of SN1987A provides strong evidence 
that the ejecta in macroscopically mixed and presumably clumpy (McCray 1993). However, 
the fraction of this escaping energy is small, and throughout this paper we will assume that 
all the radioactive energy is deposited uniformly in the ejecta. The rate of energy deposition 
by thermalized 56 Co 7-rays in the ejecta of SN1987A is given by Woosley et al. (1989): 

L 7 = 9.54 x 10 41 x exp(-t/r 56 ) x (1 - exp[-r (t/t )" 2 ]), (8) 

where L 1 is given in erg s _1 , T56 = 111.26 d is the e-folding time of 56 Co decay, and r 
is effective 7— ray optical depth of the ejecta at some fiducial time to- L 1 scales linearly 
with the mass of 56 Co in the ejecta. Therefore, the destruction rate by Compton electrons 
for species % in s _1 per particle, ko, can be written in terms of M 56 (1987A), the mass of 
56 Co produced in SN1987A, as (Liu & Dalgarno 1995, CL08): 



kcii) 



5.95 x 10 

WiXNi 



r>:; 



Mi 



r,fi 



M 56 (SN1987A) 



x exp(-t/r 56 ) x (1 - exp[-T (t/i ) 2 ]), (9) 



where iVj is the total number of particles of species i, M 56 is the mass of 56 Co in the ejecta 
of the SN being studied, and Wi is the mean energy (in eV) per ion-pair, dissociation or 
excitation for species i. Wi is defined as the ratio of the energy of the incident electron 
divided by the number of ionization, dissociation or excitation produced by collision with 
the incident electron until it comes to rest (Liu & Dalgarno 1994, Dalgarno et al. 1999). 

The effective 7— ray optical depth r(t) at time t can be written as: 

T{t) = K 5G X<j)(t) (10) 

3 Mfle 



«56 X p(t) R(t) = K 5e X 



4 7TR(t) 2 



where is the average 7-ray mass absorption coefficient in cm 2 g _1 , 0(t) is the mass- 
column density of the ejecta in g cm -2 , p[t) and R(t) are, respectively, is the mass density 
and radius of the ejecta at time t, and M# e is the mass of the helium core. The average 
mass absorption coefficient depends on the ejecta composition and the distribution of the 
radioactive material within the ejecta. For slabs of material consisting of pure He, C, O, Mg, 
Si, or Fe, k(E^ = 1.25 MeV)~ 0.056 cm 2 g _1 . For the mixed distribution of 56 Co in model 
10HMM of SN1987A, Woosley et al. (1989) derive an effective value of k 56 = 0.033 cm 2 g" 1 , 
which is the value that we will adopt for all models in this paper. Table 1 lists the value 
of t = r(t ), for t = 100 d, and the relevant parameters used in its calculation, for the 
different primordial SNe used in this study. 
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The interaction of the Compton electrons with the molecules leads to their dissociation, 
ionization and fragmentation into ionic products. The branching ratios for the different 
processes depends on Wi, the mean energies per ion-pair for a given species. Available 
values of Wi for molecules that can form in the ejecta are listed in Table 3. When data are 
not available, we just assume values similar to those for CO. Using the 56 Co mass listed in 
Table 1 for each progenitor, we calculate the time dependent rates from eq. (Q. However, 
the rate values need to be expressed under a Arrhenius temperature-dependent form as given 
by eq. We thus compose the reverse of the time-dependent temperature functions given 
by equations ([I]) and (J2J) for our modeled ejecta with our rate function given by eq. (Q, 
and derive Compton electron destruction rates as a function of ejecta temperature. We then 
fit those rates with a Arrhenius function as given by eq. ([7]). The corresponding Arrhenius 
parameters are listed in Table 3. 



3.1.2. Destruction by UV Radiation 

As they slow down in the ejecta, the high-energy Compton electrons deposit their energy 
in three channels: heating, excitation, and ionization. The ionization of the ejecta by the 
primary electrons creates secondary fast electrons, which also deposit their energy in these 
channels. The fraction of the energy going into each channel depends primarily on the 
ejecta composition, and on x e , the electron fraction of the ejecta. Kozma & Fransson (1992; 
hereafter KF92) calculated these quantities for different ejecta composition as a function of 
x e (Figs 3-5 in their paper). Of particular interest to our study is the fraction a of energy 
that is deposited in the ejecta that emerges as UV radiation. KF92 derived the value of a by 
calculating the evolution of x e for different layer compositions thus determining the fractional 
energy deposited by the high-energy electrons in the different channels as a function of time. 
They then calculated the amount of UV emission that is released by excitation and ionization 
in the different composition zones as a function of time. They found that the fraction of the 
energy that is deposited in the ejecta that emerges as UV photons with wavelengths below 
3646 A is slowly rising from a value of ~ 0.25 on day 200 to ~ 0.4 on day 1000 (Fig 8 in 
their paper). The destruction rate of molecules by UV photons (in s _1 molecule -1 ) is thus 
given by: 

x exp(-t/r 56 ) x (1 - exp[-T (t/t r 2 ]), (11) 

where a is the fraction of deposited energy re-emitted as UV photons, and Euv is the energy 
of the UV photons. According to FK92, non-thermal excitations and recombinations in the 
oxygen-rich zone result in emitting OI 1302 A and 1356 A photons. We therefore calculate 
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Euv as being the energy of a fiducial photon of wavelength 1302 A. Other terms are as in 
equation (Q. 



4. 



RESULTS 



In the following sections, we present the abundances and the mass ejected at day 1000 
of molecules produced in fully-mixed and unmixed ejecta of supernovae, excluding those that 
are dust molecular precursors. Results for gas-phase dust precursors and small clusters will 
be presented in a forthcoming paper (Cherchneff & Dwek in preparation). 



Previous studies of molecular formation in SN1987A assumed that the chemistry was at 
steady state, implying that species number densities did not vary with time (Lepp et al. 1990, 
Liu & Dalgarno 1994, 1995, Gearhart et al. 1999, Clayton et al. 2001). Such assumption 
was also made by Todini & Ferrrara (2001) and SFS04 in their studies of dust formation in 
primeval supernovae. When a small number of chemical reactions is considered, the steady 
state approximation may be valid, allowing for the direct derivation of analytical solutions 
for chemical abundances. This assumption requires that at early times, the rates of the 
chemical reactions be fast compared to the rate of density and temperature changes in the 
ejecta, so that chemical abundances would quickly reach their equilibrium value. However, 
fast chemistry does not 'a priori' ensure the validity of the steady state approximation in 
larger chemical systems. Indeed, such large systems bring a greater complexity in terms of 
chemical rates. Certain species may reach steady state abundances while others are still 
being formed or destroyed under non-equilibrium conditions. 

Figure 2 illustrates this point for the 170 M progenitor fully-mixed ejecta, by com- 
paring select reaction rates to kd yn , the inverse dynamical timescale, which is defined as: 



where v is the ejecta velocity given by eq. ([5]) and listed in Table 1, and R(t) is the ejecta 
position at a given time t after explosion. The rates depicted in the figure correspond to 
those of the major chemical processes involved in the formation and destruction of carbon 
monoxide, CO. The formation processes are the radiative association reaction (RA4) 
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and the bimolecular process (NN56) 



C + 2 



co + o. 
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The major CO destruction processes are the bimolecular process (NN70) 



Si + CO -> SiO + C, 
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the reaction with helium ions (IM13) 



He+ + CO -> C+ + O + He, 
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and the collision with 7-rays-induced Compton electrons (CED31) given by 



CO + e, 



Compton 



C + O + e, 



Compton " 
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We see from Figure 2 that prior to ~ 480 days post-explosion, all reaction rates are larger 
than the dynamical rate kd yn - A fast chemistry takes place but reactions proceed more or 
less efficiently, resulting in a non-equilibrium chemistry and a drive towards formation of 
molecules. At times greater than 480 days, the RA rate becomes smaller than kd yn and the 
reaction 'freezes' out when bimolecular processes are still active in building up molecules. 
Compton electron and He + reactions are always important in destroying CO up to 900 days 
and 700 days post-explosion, respectively. Therefore, the relevant processes to CO chemistry 
are activated and switched off at different times in the ejecta, and the chemistry is by no 
means at steady-state over the time period of interest. Consequently, molecular abundances 
show strong variations with time regardless of the initial composition and mixing in the 
ejecta. This point will be discussed in more detail below. 



We first consider the chemistry of the fully-microscopically-mixed helium cores for the 
primordial SN surrogates under study, exploring the impact of hydrogen mixing on the ejecta 
chemistry and the dependence of results on the progenitor mass. No UV destruction rates 
have been included in these models. 



4.2. Fully-mixed Ejecta 



4-2.1. Impact of Hydrogen Mixing 



Molecular abundances for the 170 M Q surrogate are shown in Figure 3 considering two 
cases: (a) when no hydrogen is mixed into the He core, and (b) when 10% of the hydrogen 
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mass of the progenitor's envelope is microscopically-mixed uniformly throughout the He core. 
In the absence of hydrogen, there exist two phases of molecular formation as illustrated in 
Figure 3a. The first phase arises at ~ 350 days when three molecules dominate the ejecta: 
silicon monoxide, SiO, silicon sulfide, SiS, and carbon monoxide, CO. For SiS and CO, the 
dominant building process is radiative association accounting for more than 94 % of the 
total formation rate, whereas for SiO, RA accounts for only ~ 56% of the formation rate. 
The extra formation pathway is the Si reaction with CO leading to the buildup of SiO from 
CO. Destruction of SiS and SiO occurs through reactions with He + at the same time that 
SiO is being formed by the destruction of CO. At times close to 440 days, silicon-based 
dust precursors (i.e., small (SiO)„ and (Si0 2 ) n clusters) start forming due to the cooler gas 
temperatures, and thus deplete Si-bearing molecules and SiO. Reaction with He + remains an 
important destruction channel for SiS, SiO and CO but for the latter, fast Compton electrons 
also contribute to the destruction. The second phase of molecular formation starts at ~ 700 
days when most of SiO is depleted into dust precursors. CO molecules keep forming from 
RA and NN reactions involving atomic oxygen and carbon. The formation of O-bearing 
species like 2 and SO are mainly triggered by NN processes involving atomic oxygen and 
sulphur. Molecular destruction is driven by NN processes as well since He + has recombined 
to neutral and the level of ionization has decreased by a factor of 100 in the ejecta. To 
summarize, the first phase of molecular formation at high temperatures is triggered by RA 
processes involving atomic species, whereas the second low-temperature phase is driven by 
NN reactions with small activation energy barriers. 

When hydrogen penetration from the progenitor envelope is allowed to a level of 10 % of 
the progenitor hydrogen mass, an active and complex NN chemistry induced by radicals such 
as OH and CH, comes to play. The radical reservoir is fed by the products of reaction between 
heavy elements and molecular hydrogen. These radicals boost at early post-explosion times 
the formation of molecules like SO, NO, and O2, the latter contributing to CO formation 
via reaction with atomic carbon. Although some of the processes mentioned in the previous 
section (i.e., reaction with He + , RA processes) are still active, the whole kinetics is dominated 
by NN processes with and without activation energy barriers. Indeed, we see form Figure 3b 
that molecular formation is already fully developed at day 300, converting a large fraction of 
the ejecta into molecules. After 700 days, CO is converted to CO2 from its reaction with OH 
and the ejecta is composed of O2, SO, CO2 and H2. The molecular masses for the dominant 
species ejected at day 1000 are listed in Table 6 for the fully-mixed ejecta as well as the 
efficiency at forming molecules defined as the ratio of the molecular mass to the helium core 
mass. For the H-free fully-mixed ejecta of the 170 M & progenitor, the dominant ejected 
species are CO and SO and the total molecular component of the ejecta equals ~ 11.3 M Q . 
Defining the molecule formation efficiency as the fraction of the ejected He-core mass that 
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is converted to molecules, we get a formation efficiency of ~ 13.7 %. As expected, its H-rich 
counterpart produces a much larger molecular mass. It chiefly forms O2, SO and CO2 and 
the total molecular mass is ~ 47 M , corresponding to a formation efficiency of ~ 57 %. 

Models assuming fully microscopically mixed ejecta are highly unrealistic, since obser- 
vations of young SN remnants show that their ejecta are unmixed. Our purpose in admixing 
H is to illustrate the paramount effect of hydrogen on molecular synthesis in the ejecta. 
Hence, any hydrogen that microscopically mixes with nearby heavy elements at the interface 
of finger structures and inhomogeneities should induce efficient synthesis of H-bearing radi- 
cals and molecules. Conversely, He + is detrimental to molecular formation, as first stressed 
by Lepp et al. (1990) in their attempt to modeling CO in SN1987A. Indeed, its attack rep- 
resents one of the principal destruction channels to molecules whereas attack by Compton 
electron always plays a minor role. This may not be the case for the unmixed ejecta, as 
discussed in the next section. 



4-2.2. Impact of Progenitor Mass 

We now turn to study the impact of the progenitor mass on the ejecta chemistry. We 
consider a massive progenitor surrogate of mass 270 M & so that we can directly compare 
results to those obtained in the previous section for the 170 M PISN progenitor. We 
also consider a CCSN with a 20 M Q progenitor to see if lower mass ejecta can efficiently 
form a molecular phase or not. As shown in Tables 1 and 4, different masses for Pop. Ill 
stars imply different initial chemical compositions, explosion energies, ejection velocities, gas 
temperatures, and number densities. The gas column depths and opacities will change too, 
resulting in the Compton electron destruction rates listed in Table 4. In both cases, no 
hydrogen penetration into the helium core is included. Molecular abundances with respect 
to total gas number density are presented in Figure 4. 

From the comparison of Figures 4 and 3 we see that the same species form in the 170 
M and the 270 M & ejecta. However, molecular formation is delayed for the more massive 
ejecta to day 550 primarily due to the higher gas temperatures. The chemical processes at 
play are identical to those mentioned in Section 4.2.1. SiO is too depleted on day 500 at 
the onset of dust precursor nucleation. As to molecular abundances, they are globally lower 
for the 270 M Q progenitor due to a less favourable initial chemical composition of the ejecta 
gas at 100 days post-explosion. Indeed, we see from Table 4 that the mass of helium (and 
thus He + ) relative to the total progenitor mass is larger for the 270 M Q surrogate compared 
to the 170 M Q progenitor. The oxygen and carbon contents are also lower, resulting in 
smaller CO abundances in the ejecta. These conditions conspire to delay the formation of 
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molecules, which consequently takes place at lower gas densities resulting in lower formation 
efficiencies. The results are summarized in Table 6 which compares the molecular yield of 
the 170 Mq and the 270 M & PISNe without and with hydrogen mixing. Therefore, low-mass 
PISNe coming from Pop. Ill progenitors should form more molecules in their ejecta than 
their very massive counterparts. 

It is of interest to compare the present results for the 170 M and 270 M surrogates 
to the study of dust formation in the fully-mixed ejecta of PISN by SFS04. Their modeling 
of dust formation accounts for the formation of CO and SiO following the formalism of 
Todini & Ferrara (2001). However, both studies derive CO and SiO masses assuming those 
molecules form at steady state, and consider only one formation channel (RA reaction) and 
one destruction channel via collision with Compton electrons. SFS04 derive a total CO mass 
of ~ 6 and 0.01 M in the ejecta of the 170 M and 260 M PISN, respectively, but they 
do not show the evolution of molecular mass with time. In our model, the CO mass ejected 
at 1000 days after the explosion of the 170 M progenitor is 5.8 M , seemingly in good 
agreement with the SFS04 value. However, we ascribe this result to pure coincidence, since 
the chemical approaches and physical models used in our studies are totally different from 
theirs. In particular, SFS04 assume much lower temperatures in their PISN ejecta than ours. 
Furthermore, their initial carbon yield is that of Heger & Woosley (2002) which is twice as 
large as our initial carbon mass taken from Umeda & Nomoto (2002), suggesting a lower CO 
formation efficiency in their model. There is also a great discrepancy between the CO yield 
calculated by SFS04 for their 260 M progenitor, and our CO yield for the 270 M progenitor, 
which is about 300 times larger. We see from Figure 4 that CO formation proceeds at times 
greater than 500 days via non-steady state chemistry. NN bimolecular processes such as 
the reaction of atomic carbon with molecular oxygen commands CO formation whereas the 
dominant destruction channel is He + attack. These mechanisms are not considered in the 
SFS04 study. All shortcomings in the SFS04 models point out the importance of using a 
kinetic approach to both processes, the synthesis of molecules and the nucleation of dust 
in the ejecta. The two are intrinsically linked, since the formation of molecules (that are 
not dust precursors) depletes the ejecta from refractory elements that would otherwise be 
included in the dust formation process. 

In the 20 M case, molecular formation patterns are quite different. On one hand, the 
lower ejecta temperatures foster molecular synthesis. On the other hand, the large helium 
content of the ejecta and the resulting higher He + abundances efficiently inhibit molecule 
production. The initial elemental composition also implies less refractory elements like Si 
or Mg available for the formation of molecules and the nucleation of metal oxides. The 
combination of the paucity of refractory elements and the elevated He + abundances at times 
less than 700 days explain the low molecular abundances, as illustrated in Figure 4. Once 
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He + has recombined at late times, the formation of molecules like CO, SO and O2 can 
proceed. We see from Table 6 that the fully-mixed ejecta of CCSNe are as efficient as their 
massive counterparts at synthesizing molecules at late post-explosion times. 

4.3. Unmixed Ejecta 

As stressed in the previous Section, it is unlikely that SN ejecta are fully microscopically- 
mixed. We thus consider totally unmixed ejecta in which each mass zone of the helium core 
has retained the stratified pre-explosive stellar composition, except for the inner most mass 
zone whose composition reflects the explosion nucleosynthesis. We further assume that the 
elements are microscopically-mixed within each mass zone. Results for the unmixed ejecta of 
two 'primordial' SNe, one 170 M® PISN and a 20 M Q CCSN, are presented in the following 
sections. The zoning and initial elemental compositions are those of Table 5, and we explore 
the impact on the chemistry of a secondary ultraviolet (UV) field as defined in Section 3.1.2. 

4.3.1. 170 M Q Ejecta: Impact of UV Radiation 

Molecular abundances normalized to the total gas number density are shown in Figure 5 
for zones 1 and 2, in Figure 6 for zones 3 and 4, and in Figure 7 for zone 5. In zone 1, the 
formation of S2 is coupled to that of SiS and FeS. The major formation process for SiS is 
the NN72 reaction S 2 + Si — ► SiS + S , whilst the RA reaction between atomic S and Si 
contributes to a lesser extent. The reverse reaction of the NN channel is actually the major 
formation process for S2. However, since the initial mass of Si is four times larger than that 
of S, the net formation of SiS is always more efficient than that of S2. After ~ 350 days, 
the reservoir of S2 is depleted through the continuous formation of SiS and the simultaneous 
formation of FeS, the molecular precursor to iron sulfide clusters, through the reaction of Fe 
with S2 (see Figure 5a). As a result, SiS is the most abundant molecules formed in zone 1. 

In zone 2, most molecules are rapidly formed from RA processes whereas thermal frag- 
mentation is the dominant destruction process at early times, due to the high gas tempera- 
tures. O2 and SiO formation are coupled, for O2 is partly destroyed by the NN69 reaction 
with atomic Si to form SiO. At ~ 250 days post-explosion, SiO abundance decreases due 
to nucleation of silica dust precursors such as (SiO) n , resulting in a sharp increase of O2 
abundance. As SO and CO are also formed in the NN reaction of O2 with S and C atoms 
(reactions NN76 and NN56, respectively), their abundances show a steep rise with increasing 
O2. The low amount of CO formed in the gas is due to the initial low carbon content of this 
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mass zone. 

The chemical processes pertaining to zone 2 also apply to zone 3 and the resulting 
molecules are similar to those of zone 2, as seen in Figure 6a. The differences in abundance 
variation and magnitude with time are due to the different initial chemical composition of 
zone 3. Specifically, the higher carbon initial mass yield drives the formation of CO via RA 
process at early times, depleting some of the atomic oxygen available to the formation of O2 
and resulting in higher CO abundances and lower O2 and SO abundances over time than in 
zone 2. 

Despite the fact that zone 4 is oxygen-rich, it is characterized by high initial carbon mass 
compared to zones 2 and 3 (see Table 5). However, it is expected that its initial chemical 
composition precludes formation of carbon-rich molecules other than CO as its C/O ratio is 
~ 0.29. However, we see from inspection of Figure 6b that very early on, C2 and C3 form in 
large amounts. This carbon chain formation is a consequence of the thermal fragmentation 
of CO at very high temperatures and its rapid conversion to C2. Indeed, fragmentation 
studies of CO in shock tube experiments at temperatures between 5000 K and 18000 K 
show that C2 and electronically-excited CO* are always detected along with the collisional 
destruction of CO in the high temperature post-shock gas (Fairbairn et al. 1968, Appleton 
et al. 1970, Hanson 1973). This is readily explained by the following chemical mechanism 
(labelled TF13, NN57 and TF18 in Table 9): 

CO + M C + + M, (18) 

C + CO <-> C 2 + O, (19) 

and 

C 2 + M^C + C + M. (20) 

where M is the colliding buffer gas (chosen as Ar, O, C, or CO in shock tube experiments), 
and all reactions and their reverse processes occur simultaneously. We see that along with the 
destruction of CO by collision with M, its rapid conversion to C2 occurs through reaction with 
atomic carbon [eq. fll9p ]. Chemical rates characterizing these reactions have high activation 
energy barriers and thus only proceed at high temperatures. Indeed, at temperatures less 
than 5000 K, rates for reactions (18), (19) and (20) become small. As to C3, the build-up 
of the end-product of the carbon chains is due to its formation from reaction of two C2 
molecules which rapidly depletes C2 chains from the gas, as seen in Figure 6b. A carbon 
dust mass yield can not be derived from the C3 abundance profile, for the formation of 
C-bearing large chains and rings has not yet been included in the chemical network. It is 
expected that C3 will be transformed in larger carbon chains and that atomic oxygen in 
the gas will partly destroy those chains through the formation of CO. Therefore, carbon 
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chain abundances should be smaller than values for C3 derived in this study. However, the 
quick conversion of CO into C2 will keep triggering the build-up of carbon dust molecular 
precursors in this O-rich mass zone. Two rates are reported in the literature for the C + CO 
reaction: a fast rate calculated by Fairbairn et al. (1968) and Hanson et al. (1973), which 
we use in our ejecta model, and a value listed in the UDFA06 database but with no traceable 
reference (Woodhall et al. 2007). The latter is ~ 10 times smaller than the former. We 
ran the full model for zone 4 including the lowest rate and find that the C3 abundance at 
day 1000 is decreased by a factor of 3, the CO abundance is increased by a factor of 6, and 
the O2 abundance is decreased by 30% compared to the values illustrated in Figure 6b. In 
both cases however, the quick CO conversation to C2 via the C + CO reaction observed 
in thermal fragmentation experiments is a 'natural' mechanism to produce carbon-bearing 
species and carbon dust in a hot, oxygen- rich environment where carbon is initially in atomic 
form. When CO conversion to C2 is not considered, no C2 forms in large enough quantities 
to provide significant amounts of the end-product carbon chain C3. 

To account for the formation of carbon dust in CCSNe, Clayton et al. (1999, 2001) 
proposed that CO dissociates in a fully-mixed ejecta from collisions with Compton electrons 
and reactions with He + . The chemistry is assumed to be at steady state and CO dissociation 
creates a pool of carbon atoms from which carbon chains and solid clusters can nucleate 
and condense. However, we see from Section 4.1 that under non-steady state conditions, 
the two mechanisms invoked by Clayton et al. do not succeed in destroying enough CO 
to provide free carbon atoms. Indeed, for the three hydrogen-free mixed ejecta, CO is 
always the dominant species synthesized in the gas. Furthermore, He + being a predator 
to molecules at early times, CO formation is delayed to post-explosion times characterized 
by temperatures low enough to preclude the quick CO conversion to C2 from operating. 
We conclude that fully-mixed SN ejecta can not form free atomic carbon and carbon-based 
solids in general. Conversely, unmixed SN ejecta, in particular He-free zones where oxygen 
is more abundant than carbon, can generate carbon chains from the quick conversion of 
CO through the C + CO reaction operating at high temperatures. As to prevalent chemical 
formation processes in zone 4, we find that both RA reactions and NN bimolecular routes are 
active at high temperatures to forming species which are destroyed by thermal fragmentation 
reactions. Molecular formation and destruction in the intermediate temperature regime (1000 
K - 5000 K) are too governed by NN bimolecular reactions when destruction via collisions 
with Compton electron also participate but to a much less extent. 

As to zone 5, we see from Figure 7 that the major molecule to form is carbon monoxide. 
In the outer zone, the formation of He + at early times triggers the simultaneous formation 
of ions like + and C + . NN processes are then not as important as they are in other zones 
for the formation of molecules. This fact, combined with the destruction of species by He + , 
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hampers the effective formation of molecules over most of the post-explosion times. CO, and 
to a less extent O2, thus form at late times after He + recombination. 

The molecular masses ejected at day 1000 are summarized for all zones in Table 7 and 
illustrated in Figure 8. The efficiencies listed are defined as the ratio of the molecular mass 
formed per zone to the zone mass. The value in the last column corresponds to the total 
molecular mass formed in the unmixed ejecta divided by the summed mass of the zones. 
In the unmixed case, molecules reflect the chemical composition of the zone in which they 
form. Dioxygen O2 is by far the most abundant species in the ejecta and form efficiently in 
three mass zones when SiS too forms in large amount but is confined to the innermost region 
of the He core. The final molecular content of the ejecta equals ~ 37.1 M Q , and the total 
formation efficiency is 45.07 %, a very high value. It means that massive PISNe are effective 
at converting almost half of their initial atomic content into a molecular phase. This value 
is much larger than that for the fully-microscopically mixed case without hydrogen mixing. 
This result once again illustrates the primary role of He + in impeding the formation of 
molecules. 

The impact of the UV field as described in Section 3.1.2 is explored for zone 4 which 
produces the largest amount of molecules. The results on molecule abundances are shown in 
Figure 9. Formation and destruction processes for the dominant species are similar to those 
when no UV field is considered. While thermal fragmentation and NN bimolecular reactions 
remain the dominant destruction processes at high temperatures, molecular destruction by 
UV becomes effective at times greater than ~ day 350, reducing the abundances of some 
important molecules. Specifically, O2 is affected as its formation is postponed to later times, 
as seen in Figure 9. The final O2 abundance at day 1000 is decreased by 30 % compared to 
its value when no UV radiation is considered. However, the overall major chemical processes 
involved in molecular synthesis remain identical. We also consider UV photodissociation 
in our fully-mixed SN ejecta and find that the impact of UV photodissociation is minor 
compared to the destruction of chemical species by He + . We thus conclude that UV photo- 
processes if present have some minor impact as destruction channels to molecules and that 
their effect on the overall chemistry of the ejecta is limited. 



4.3.2. 20 M & Ejecta 

We now turn to studying the unmixed ejecta of a zero-metallicity, 20 M Q progenitor 
core-collapse supernova. Results for zones 1, 2 and 3 as defined in Table 5 are presented in 
Figures 10 and 11. These zones provide most of the molecules in the ejecta as seen from 
Table 8, which gives the total molecular masses formed in each zone at day 1000. These 
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mass yields are illustrated in Figure 12. 

The chemical processes which occur in zone 1 are similar to those at play in zone 1 of 
the 170 M Q unmixed ejecta, resulting in the ejection of large amounts of SiS a day 1000. 
However, large amounts of oxygen and carbon compared to silicon and sulphur are initially 
presented in the initial composition of the 20 M CCSN, resulting in an active carbon 
chemistry. For example, CS forms in large quantities from the reaction of atomic C with S2 
and CO abundances are too enhanced compared to zone 1 of the 170 M & ejecta. 

The chemistry of zone 2 is similar to that at play in zone 2 of the 170 M PISN, that is, 
early formation of most molecules by RA reactions and destruction by thermal fragmentation, 
and coupled chemical processes for the formation of O2, SiO, CO and SO at later times. The 
dominant species ejected at day 1000 are O2, SO and CO, as illustrated in Figure 10b. 

Zone 3 can be compared to zone 4 of the 170 M & progenitor both zones are 

characterized by a C/O ratio less than 1 (0.33 and 0.29 for the 20 M & zone 3 and the 
170 M e zone 4, respectively) and are helium-free. Similar chemical processes are effective at 
building up molecules in both zones. The rapid conversion of CO to C2 once again triggers 
the formation of carbon chains in zone 3 but the overall process is less efficient than for 
the 170 Mq progenitor owing to the lower temperatures in the ejecta. This is illustrated in 
Figure 10b where CO always remains more abundant than the carbon-chain end product C3, 
and is gradually converted into 2 at late post-explosion times by its reaction with atomic 
oxygen. 

In zone 4, characterized by a large He mass and a C/O ratio of 29.5, the mass of 
molecules formed is negligible as seen from Table 8. This is primarily due to the initial 
chemical composition of the zone where more than 98% of the mass is helium while C and 
O only represent 1.7% and 0.08% of the zone mass, respectively. Molecules once formed are 
chiefly destroyed by He + at any times in the ejecta owing to its overwhelming presence. 

Table 8 and Figure 12 show that low-mass, zero-metallicity progenitors are almost as ef- 
ficient at forming a molecular phase in their ejecta than their massive counterparts. The total 
efficiency at forming molecules is 35.65%, a slightly lower value than that of the 170 M case. 
This lower efficiency is primarily due to the 20 M Q zone 4, which, as discussed above, does 
not form molecules. 
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5. SUMMARY AND DISCUSSION 



We have investigated the chemistry of the ejecta of Pop. Ill progenitor SNe and find 
that copious amounts of molecules form in these inhospitable environments. Of particular 
importance are the following points: 



• As already stated by CL08, the chemistry in SN ejecta is not at steady state from 100 to 
1000 days after explosion. New chemical channels involving neutral-neutral processes 
prevail over radiative association reactions, Compton electron destruction routes and 
photodissociation by ambient ultraviolet photons. Ion-molecule reactions play a role 
at late times when the ejecta is cool and diffuse. Our results for relevant molecules 
like CO or SiO disagree with existing studies due to the fact that the chemistry is not 
at steady state. We find that the chemistry of Pop. Ill SN ejecta, owing to the large 
ranges of temperatures and densities spanned over relatively short times, is complex, 
manifold, and conducive to molecule synthesis. 

• A new pathway to the formation of carbon chains is active in the oxygen-rich mass zone 
of the unmixed ejecta and is identified as the CO conversion to C2 via collision with 
C. This fast conversion is usually observed in the thermal fragmentation of carbon 
monoxide in high temperature shock tube experiments. When this conversation is 
suppressed, no carbon-bearing molecules and chains are formed. Thus, this conversion 
channel triggers the formation of carbon chains and dust in an oxygen-rich gas. It is 
then relevant to any gaseous O-rich environment characterized by high temperatures 
(T > 5000 K) and a large atomic carbon fraction. 

• The present results are extremely sensitive to mixing in the ejecta. We find that the 
injection of hydrogen from the progenitor envelope in fully-mixed ejecta boosts molec- 
ular synthesis via the formation of radicals like OH. This result too applies to unmixed 
ejecta. On the other hand, we show that helium severely hampers the formation of 
molecules through He + attack. Therefore, the detection of molecules in SN ejecta 
brings evidence for the non-mixing of helium with other elements in the ejecta gas. 

• The results of our calculations show that the fully-mixed 170 M & and 270 M & pro- 
genitors produce 11 M and 3.2 M Q of molecules, respectively. Therefore, a larger 
progenitor mass does not imply a larger molecular content of the ejecta. This is chiefly 
due to the harsh physical conditions encountered in the ejecta of the 270 M & progen- 
itor and to its initial chemical composition. Indeed, although the 270 M model is 
characterized by larger heavy element masses with respect to its 170 M counterpart, 
its helium mass is almost three times greater than that of the 170 M & case, implying 
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efficient destruction of molecules. The admixing of 10% of the hydrogen present in 
the 170 M & progenitor envelope into the fully-mixed ejecta dramatically increases its 
molecular yield to ~ 47 M Q . The more realistic unmixed ejecta of a 170 M & progeni- 
tor supernova synthesizes ~ 37 M Q of molecules at post-explosion day 1000, which is 
significant. About half of the initial elemental content of the He core is converted into 
molecules. The most abundant species by mass is O2 followed by SiS, CO and SO. 
Its 20 M counterpart produces ~ 1.2 M . 2 is the dominant species followed by 
CO, SiS, and SO. This lower efficiency at forming chemical species for the low mass 
CCSN is due to the existence of its extended helium-rich outer zone in which molecular 
synthesis is suppressed. 

• The present results hold for Pop. Ill SNe but the large mass yields of molecules 
formed in their ejecta address the possibility of potential observational detection of 
new molecules in nearby SN ejecta. As stated above, our primordial 20 M Q progenitor 
forms O2, CO, SiS, and SO shortly after explosion. CO and SiO have already been 
detected at IR wavelengths in several SN ejecta. However, molecules like SiS and O2 
are tracers of stratified, unmixed ejecta while microscopic mixing with hydrogen pro- 
duces tracer species like CO2, OH and H2O. Search for those chemical species should 
be undertaken at IR and submillimetre wavelengths. 

The large amounts of molecules synthesized in the PISN ejecta are exposed to a harsh 
environment generated by the PISN blast wave. The shock expanding into the ambient 
circumstellar /interstellar medium generates fast particles that are accelerated to cosmic rays 
energies. These penetrate the cavity generated by the blast wave subjecting the ejecta to the 
constant bombardment by high energy electrons and nuclei. The pressure of the material 
that is shock-heated by the advancing blast wave will generate a reverse shock that will 
move into the expanding PISN ejecta (e.g., Nozawa et al. 2007). Locally, observations of 
young supernova, such as Cas A, show that their ejecta is very clumpy, the clumps consisting 
of X-ray filaments with densities ~ 1 — 10 cm -3 , and optical- and IR-line emitting knots 
with densities of ~ 10 3 — 10 4 cm -3 (Fesen et al. 2006, Smith et al. 2009). Therefore, 
the fate of the ejected molecules will depend on their environment. Low density filaments 
are heated by the reverse shock to X-ray emitting temperatures. Chemical species are not 
likely to survive this harsh environment. The dust and molecules inside the optical- and 
IR-line emitting clumps will encounter a much slower shock. Their fate will depend on 
the relative timescales of many different processes operating in the cavity of the young 
remnant: the radiative cooling time of the shocked clump, the timescales for the heating 
of the clump by thermal conduction and ambient cosmic rays, the evaporation timescale 
of the clump, and the timescale for the development of various instabilities that can lead 
to its fragmentation and subsequent evaporation. Any surviving molecules will be further 
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subjected to the general diffuse interstellar UV radiation field generated by the Pop III stars. 
It is therefore very unlikely that the molecules synthesized in the Pop III SN ejecta will have 
any global cosmo logical impact. However, they can have a significant local impact. The 
expanding PISN blast wave will generate during the radiative phase of its evolution a cold 
dense shell. This shell may be subject to various instabilities that can cause its collapse, 
forming the next generation of stellar objects (MacKey et al. 2003, Schneider et al. 2006). 
These stars, commonly referred to as Pop 2.5 stars, will form out of a gas that contains an 
admixture of the heavy elements, molecules, and dust that formed in the PISN ejecta. The 
mass of these stellar objects will depend on their ability to fragment into smaller structures, 
which is greatly facilitated by the cooling rate of the gas via atomic and molecular processes, 
and by the conversion of the cloud's internal energy to infrared emission by dust (Bromm & 
Larson 2004). The survival of the molecules and dust, and their effect on the formation of 
stars in this propagating star formation scenario will be explored in a subsequent paper. 

IC would like to acknowledge support from the Swiss National Science Foundation 
through a Maria Heim-Vogtlin fellowship. 
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Table 1. Primordial supernova parameters used in this study (adapted from UN02, Heger 

& Woosley 2002, and NK03.) 





20 Mq 


170 Mq 


270 Mq 


Eexplosion(Ergs) 


lxlO 51 


2xl0 52 


8xl0 52 


He core mass (Mq) 


5.8 


82.3 


129 


M( 56 Co) (Mq) 


0.07 


3.6 


9.8 


v (Km s~ 1 ) 


2242 


3439 


5458 


To (K) 


18000 


21000 


31500 


t (g cm" 2 ) 


23.80 


146.07 


90.89 



Table 2. Species included in the ejecta chemical models 



Atoms 


II 


He 





c 


Si 


S 


Mg 


Fe 


Al 






Diatomic species 


Hi- 


OH 


o 2 


CO 


SiO 


SO 


NO 


MgO 


FeO 


AlO 


C 2 




CS 


CN 


SiH 


SiC 


Si 2 


SiS 


SiN 


SH 


N 2 


NH 


MgS 




Fe 2 






















Tri-atomic species 


H 2 


H 2 S 


HCN 


CH 2 


C 2 H 


HCO 


c 3 


co 2 


ocs 


OCN 


SiC 2 




Si 3 


S1O2 


so 2 


N0 2 


Fe 3 














4-atom species 


Si 2 2 


Mg 2 2 


Mg 2 S 2 


Fe 2 2 


Fe 2 S 2 


H 2 CC 


Fe 4 


Si4 


C 3 H 


CH 3 




> 5-atom species 


Si 3 3 


Si 2 4 


Si 3 6 


C 3 H 3 


C4H4 


C 6 H 5 


C 6 H 6 










Ions 


H+ 


H _ 


Hc+ 


0+ 


Si+ 


S+ 


Mg+ 


Fe+ 


A1+ 


H* 


H 3 + 




HeH+ 


ct 


CO+ 


SiO+ 


SO+ 


H 2 0+ 


HCO+ 
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Table 3. Compton electron-induced reactions, corresponding mean energy per ion pair Wi 
and Arrhenius coefficient A as a function of ejecta model. 



Species 


Reaction 


W, (eV) 


A - 20 Mq a 


A - 170 M Q a 


A - 270 Mq a 


Reference 


CO 


-> o+ + C 


768 


1.1610X10" 7 


9.4576xl0~ 7 


1.6741 xlO -6 


Liu & Dalgarno (1995) 




-> c+ + o 


247 


3.6100xl0~ 7 


2.9406X10" 6 


5.2053xl0~ 6 






-» c + o 


125 


7.1333xl0~ 7 


5.8107xl0~ 6 


1.0286X10 -5 


» 




-> CO+ + e- 


34 


2.6225X10" 6 


2.1363xl0~ 5 


3.7815X10- 5 




O 


-> 0+ + e- 


46.2 


1.9300X10" 6 


1.5722X10" 5 


2.7829xl0~ 5 




C 


-* C+ + c" 


36.4 


2.4496xl0~ 6 


1.9954xl0 -5 


3.5321xl0~ 6 




SiO 


-> 0+ + Si 


678 


1.3158X10" 7 


1.0719X10 -6 


1.8973X10" 6 






-» Si+ + o 


218 


4.0913xl0~ 7 


3.3327X10" 6 


5.8993X10- 6 






-» Si + o 


110 


8.0844xl0~ 7 


6.5855xl0~ 6 


1.1657X10 -8 


» 




-> SiO+ + e- 


30 


2.9722X10" 6 


2.4211X10" 5 


4.2857X10" 5 




N 2 


— > N+ + N 


264 


3.3812X10" 7 


2.7543xl0" e 


4.8755X10" 6 


Khare & Kumar (1977) 




-* N+ N 


133.5 


6.6813X10- 7 


5.4425 xl0~ 6 


9.6339X10- 6 






N+ + e- 


36.3 


2.4564 xlO -6 


2.0009xl0~ 5 


3.5419xl0~ 5 




H 


-> H+ + e- 


36.1 


2.4700 xlO -6 


2.012xl0~ 5 


3.5615X10" 6 


Dalgarno Yan Liu (1999) 




H* (n=2) 


26.6 


3.3521 xlO" 6 


2.7306xl0~ 5 


4.8335 XlO" 8 




He 


-> He+ + e- 


46.3 


1.9258X10 -6 


1.5688X10 -5 


2.7769xl0~ 5 




H 2 


— > H+ + H 


820 


1.0874X10" 7 


8.8578xl0~ 7 


1.5679xl0" 6 






-> H + H 


77 


1.1580X10 -6 


9.433xl0~ 6 


1.6697xl0~ 8 


?? 




-> H+ + e+ 


37.7 


2.3651X10" 6 


1.9266X10 -5 


3.4103xl0" 5 


J5 



a Arrhenius forms for kc (see text): A X exp(-2976.5/T) (20 Mq progenitor) - A X exp(-3464.1/T) (170 Mq progen- 
itor) - A x cxp(-5376/T) (270 Mq progenitor) 



Table 4. Initial (post-explosive) chemical composition in units of M© for fully-mixed primordial SN ejecta without 
hydrogen mixing. The initial mean molecular weight /io(gas) is given in g mole -1 . 





w(gas) 


He 





Si 


S 


Mg 


Fe 


C 


Al 


Ne 


Ar 


N 


20 M© 


5.63 


3.59 


1.55 


9.82 (-2) 


4.12 (-2) 


7.07 (-2) 


7 (-2) 


0.26 


4.79 (-4) 


0.12 


6.9 (-3) 


2.69 (-4) 


170 MO 


18.11 


1.96 


44.23 


16.16 


8.66 


1.94 


3.63 


2.30 


2.0 (-2) 


1.19 


1.42 


1.0 (-2) 


270 MO 


19.31 


5.5 


44.31 


26.95 


15.78 


4.78 


16.14 


1.89 


8.62 (-2) 


4.70 


2.60 


1.26 (-2) 



Table 5. Initial (post-explosive) chemical composition in units of MO for unmixed primordial SN ejecta. The initial 



mean 


molecular 


weight fio(gas) 


is given 


(in g 


mole 1 ) as 


well as 


the C/O ratio of each zone. 






Zone 


M0 (gas) 


c/o 


He 





Si 


s 


Mg 


Fe 


c 


Al 


Cr 


Co 


Ni 


Unmixed 20 M 


Zone 1 (2.4-3 M©) 


30.21 


0.013 





6 (-4) 


0.39 


0.138 


6 (-6) 


0.048 


6 (-6) 


1.2 (-7) 


3 (-3) 


4 (-4) 


2 (-4) 


Zone 2 (3-3.6 MQ) 


16.87 


0.0013 





0.52 


0.0358 


0.0072 


0.0363 





4.98 (-4) 


4.2 (-5) 











Zone 3 (3.6-4.95 MQ) 


15.03 


0.33 





1.08 


2.7 (-5) 


6.75 (-7) 


4.725 (-3) 





0.266 


4.05 (-6) 











Zone 4 (4.95-5.85 M©) 


4.05 


29.47 


0.884 


S.84 (-4) 


9.0 (-9) 





2.25 (-5) 





1.512 (-2) 














Unmixed 170 M 


Zone 1 (0-20 M0) 


29.13 


0.066 





3.5 (-5) 


13.2 


4.0 


5.3 (-5) 


0.35 


1.8 (-6) 


1.8 (-10) 


3.5 (-3) 


1.2 (-3) 


3.5 (-3) 


Zone 2 (20-40 M0) 


17.29 2.9 xlO" 5 





16.5 


2.76 


0.4 


0.32 





3.6 (-4) 


5 (-4) 











Zone 3 (40-55 M©) 


16.76 


0.03 





13.1 


0.615 


3 (-2) 


1.22 





3 (-2) 


1.5 (-2) 











Zone 4 (55-78 M0) 


15.17 


0.29 





18.6 


1.84 (-3) 


1.38 (-6) 


0.299 





4.07 


4.6 (-4) 











Zone 5 (78-82 M0) 


10.46 


0.56 


0.596 


2.4 


4 (-6) 


2.4 (-7) 


3.6 (-3) 





1.0 


4 (-8) 
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Table 6. Mass yields of most important molecules ejected at day 1000 for fully-mixed 

ejecta with and without hydrogen mixing. a,fc 



Molecules 


170 M Q 


170 Mq 


270 M 




20 M C 




no hydrogen 


10 % hydrogen 


no hydrogen 


no hydrogen 


CO 


5.61 


0.18 


3.22 




0.63 


SO 


5.61 


13.28 


1.68 xlO" 


-4 


6.35 xlO" 2 


o 2 


5.24 xlCT 2 


24.71 


2.67 xlO" 


-4 


8.79 xlO" 2 


co 2 


7.20 xl0~ 5 


8.49 


3.45 xl0~ 


-8 


1.48X10" 4 


H 2 





0.12 










N 2 


1.03 xlO" 2 


3.42 xl0~ 3 


5.70 xlO" 


-4 


2.73 xl0~ 4 


NO 


5.21 xlO" 4 


1.62 xl0~ 3 


1.88 xlO- 


-7 


2.47 xlO" 6 


OH 





1.31 xlO -3 










Total 


11.29 


46.80 


3.22 




0.78 


Efficiency 


13.71 % 


56.86 % 


2.49% 




13.42 % 



a Mass yields are in Mq 

b The efficiency is defined as the ratio of the molecular mass to the He core mass 
C A mass cut of 2.4 Mq is assumed as in Nozawa et al. (2003) 
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Table 7. Mass yields of most important molecules ejected at day 1000 for the unmixed 
ejecta of the 170 M progenitor without hydrogen mixing. a ' b 





Zone 1 


Zone 2 


Zone 3 


Zone 4 


Zone 5 


Zone 1-5 


Zone mass 


(20 M ) 


(20 M Q ) 


(15 M ) 


(23 M ) 


(4.3 M ) 


(82.3 M©) 


Zone C/O 


0.066 


2.9xl0" 5 


0.03 


0.29 


0.56 




o 2 





6.24 


5.98 


13.80 


9.28 xl0~ 5 


26.00 


SiS 


7.36 








8.79 xl0~ 6 





7.36 


CO 





8.39X10" 1 


6.98 xlO" 2 


9.95 xlO" 1 


2.03 


3.09 


so 


9.06xl0~ 5 


5.97X10" 1 


4.47 xl0~ 2 








0.64 


Total mass 


7.36 


6.84 


6.10 


14.80 


2.03 


37.09 


Efficiency 


41.81% 


34.18 % 


40.64 % 


64.22 % 


47.21 % 


45.07 % 



a Mass yields are in Mq 

b The efficiency is defined as the ratio of the molecular mass to the zone mass 



-35 - 



Table 8. Mass yields of most important species ejected at day 1000 for the unmixed ejecta 
of the 20 M Q progenitor without hydrogen mixing. a ' b 





Zone 1 


Zone 2 


Zone 3 


Zone 4 


Zones 1- 4 


Zone mass 


(0.6 M Q ) C 


( 0.6 M ) 


(1.35 M Q ) 


(0.9 M ) 


(3.45 M ) 


Zone C/O 


0.013 


0.0013 


0.33 


29.47 




2 





0.29 


0.45 





0.74 


CO 





1.15XKT 3 


0.27 


3.63 xl0~ 5 


0.27 


SiS 


0.21 











0.21 


so 


9.43X10" 4 


1.03X10" 2 


9.49 xlO" 7 





0.01 


Total mass 


0.21 


0.3 


0.72 





1.23 


Efficiency 


35.37% 


50.14% 


53.41 % 


0.004 % 


35.75 % 



a Mass yields are in M 

b The efficiency is defined as the ratio of the molecular mass to the zone mass 
C A mass cut of 2.4 M is assumed for Zone 1 as in Nozawa et al. (2003) 
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Fig. 1. — The evolution of the gas parameters as a function of time for the different SN ejecta 
(taken from Nozawa et al. 2003) and assumed to be independent of mass zone within the He 
core, a) Temperature where the 270 M Q profile has been rescaled by a factor 1.5 compared 
to that of the 170 M ejecta to account for the larger explosion energy; b) Number density. 




Fig. 2. — The time dependence of the rates of the major processes leading to the formation 
of CO (eqs. (fl3j) - (JTT|) is compared to the inverse of the dynamical time scale of the ejecta 
kdyn for the fully mixed 170 M & case without hydrogen diffusion. The figure shows when 
reactions proceeds and freeze out, leading to an active chemistry far from Steady-State. 
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Fig. 3. — The evolution of molecular abundances normalized to total gas number density for 
the fully mixed ejecta of the 170 M Q progenitor when no hydrogen mixing is considered (a), 
and when 10 % of the hydrogen mass of the progenitor envelope is microscopically mixed to 
the helium core (b). The figure illustrates the sensitivity of the chemistry to the presence of 
hydrogen in the ejecta. 
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Fig. 4. — The evolution of the molecular abundances normalized to total gas number density 
for the fully mixed ejecta: a) 270 M Q progenitor; b) 20 M progenitor. No hydrogen mixing 
is considered in both cases. 
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Fig. 5. — The evolution of the molecular abundances normalized to total gas number density 
for the 170 M Q unmixed ejecta: a) The Si/S/Fe-rich zone 1; b) The O/Si/Mg/S-rich zone 2. 
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1E+00 
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Post-explosion time (days) Post-explosion time (days) 

Fig. 6. — The evolution of the molecular abundances normalized to total gas number density 
for the 170 M unmixed ejecta: a) The O/Mg/Si-rich zone 3; b) The O/C/Mg-rich zone 4 
in which the rapid conversation of CO to C 2 take place despite a C/O ratio of 0.29. 
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Fig. 7. — The evolution of the molecular abundances normalized to total gas number density 
in the O/C/He-rich zone 5 of the 170 M unmixed ejecta. CO is the only molecule to form 
due to the destruction of molecules by He + . 
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Fig. 8. — The mass yield (in M Q ) of molecules ejected at day 1000 for the unmixed ejecta 
of the 170 M & progenitor. 
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Fig. 9. — The evolution of the molecular abundances normalized to total gas number den- 
sity for zone 4 of the 170 M Q unmixed ejecta where UV radiation is included. Molecular 
abundances are quasi similar to those in Figure 6 except for O2 whose formation is delayed 
to later times due to UV destruction. 
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Fig. 10. — The evolution of the molecular abundances normalized to total gas number density 
for the 20 M unmixed ejecta: a) The Si/S/Fe-rich zone 1; b) The O/Si/Mg/S-rich zone 2. 
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Fig. 11. — The evolution of the molecular abundances normalized to total gas number density 
for the O/C/Mg-rich zone 3 of the 20 M unmixed ejecta. The rapid conversation of CO to 
C2 takes place despite a C/O ratio of 0.33. 
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Fig. 12. — The mass yield (in M Q ) of molecules ejected at day 1000 for the unmixed ejecta 
of the 20 M progenitor. 
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Table 9. Chemical processes considered in the present study. Chemical reaction rate 
coefficients are given according to equation ([7]). The Compton electron destruction and UV 
photodissociation reactions are listed for the 170 M Q progenitor. Rate values for other 

progenitors can be derived from Table 3. 







Chemical 


processes 


A 






Reference a 


TERMOLECULAR 


3B1 


H + 


H + M 


► 


H 2 + M 


6.84xl0" 33 


-1 





NIST 


3B2 


H + 


C + M 


> 


CH + M 


l.ooxicr 33 








E 


3B3 


H + 


O + M 


> 


OH + M 


4.36X1CT 32 


-1 





NIST 


3B4 


H + 


OH + M 


> 


H 2 + M 


2.59X1CT 31 


-2 





NIST 


3B5 


H + CN + M 


► 


HCN + M 


8.63X1CT 30 


-2.2 


566.5 


NIST 


3B6 


H + 


CO + M 


> 


HCO+ M 


5.29X10" 34 





370.4 


NIST 


3B7 


H + 


C 2 + M 




C 2 H+ M 


l.ooxicr 33 








E 


3B8 


H + 


C 3 + M 




C 3 H+ M 


l.OOxlO" 33 








E 


3B9 


O + 


C + M 




CO + M 


2.14X10" 29 


-3.08 


-2114.0 


UDFA06 


3B10 


O + 


O + M 




2 + M 


9.26X1CT 34 


-1 





NIST 


3B11 


O + 


S + M 




SO+ M 


9.26X1CT 34 


-1 





E as 3B9 


3B12 


O + 


N + M 




NO+ M 


5.46X10" 33 





155.1 


NIST 


3B13 


O + 


CO + M 




C0 2 + M 


1.20xl0~ 32 





2160.0 


NIST 


3B14 


O + 


Si + M 




SiO+ M 


2.14xlCr 29 


-3.08 


-2114.0 


E as 3B9 


3B15 


C + 


C + M 




C 2 + M 


l.OOxlO" 33 








E 


3B16 


C + 


C 2 + M 




C 3 + M 


1.00xl0~ 33 








E 


3B17 


C + 


S + M 




CS+ M 


2.14X10" 29 


-3.08 


-2114.0 


E as 3B9 


3B18 


C + 


N + M 




CN+ M 


9.40X10" 33 








NIST 


3B19 


Si + 


H + M 




SiH+ M 


l.ooxicr 33 








E 


3B20 


Si + 


N + M 




SiN+ M 


9.40xl0~ 33 








E as 3B17 


3B21 


Si + 


S + M 




SiS+ M 


2.14X10" 29 


-3.08 


-2114.0 


E as 3B9 


3B22 


s + 


S + M 




S 2 + M 


2.76X10" 33 








NIST 


3B23 


N+ N + M 




N 2 + M 


1.25xKT 32 








NIST 


3B24 


CO 


+ CH + M 




HC 2 0+ M 


2.80X10" 34 


-0.4 





NIST 


3B25 


H + 


CO + H 




0+ CH 2 


l.OOxlO" 33 








E 


3B26 


H + 


C0 2 + H 




2 + CH 2 


l.OOxKT 33 








E 


3B27 


H + 


C 2 H 2 + H 




CH 2 + CH 2 


1.00xl0~ 36 








E 


3B28 


CO 


+ OH + H 




2 + CH 2 


LOOxlO -33 








E 


THERMAL FRAGMENTATION 


TF1 


H 2 - 


1- M 




H + H+M 


2.54X10" 8 


-0.1 


52555.6 


NIST 


TF2 


CH 


+ M 




C + H+M 


3.16X10" 10 





33700.0 


NIST 


TF3 


OH- 


h M 




O + H+M 


4.00xl0~ 9 


-0.1 


50000.0 


NIST 


TF4 


SiH 


+ M 




Si+ H+M 


3.16X10- 10 





33700.0 


E as TF2 


TF5 


H 2 + M 




OH + H+M 


5.80xl0~ 9 





52920.0 


NIST 


TF6 


HCN+ M 




H+CN+M 


2.08xl0~ 8 





54630.0 


NIST 


TF7 


CH 2 


+ M 




CH+H+M 


6.64xl0- 9 





41852.0 


NIST 


TF8 


CH 2 


+ M 




H 2 +C+M 


2.66X10" 10 





32230.0 


NIST 


TF9 


C 2 H+ M 




C 2 +H+M 


3.75xl0~ 10 





50040.0 


NIST 


TF10 


C 3 H+ M 




C 3 +H+M 


l.OOxlO" 10 





48600.0 


E 


TF11 


HC 2 


0+ M 




CO+CH+M 


1.08x10-* 





29585.0 


NIST 


TF12 


2 4 


■ M 




O+O+M 


5.17X10- 10 





58410.0 


NIST 


TF13 


CO- 


- M 




C+O+M 


4.40X10- 10 





98600.0 


Appleton 1970 


TF14 


SO-f 


- M 




S+O+M 


6.61X10- 10 





53885.0 


NIST 


TF15 


SiO+ M 




Si+O+M 


4.40X1Q- 10 





98600.0 


E as TF13 
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Table 9 — Continued 







Chemical processes 


A 


V 


E a 


Reference a 


TF16 


NO- 


- M 




N+O+M 


4.10X10" 9 





75380.0 


NIST 


TF17 


CO2 


+ M 




CO+O+M 


8.02X10" 11 





26900.0 


NIST 


TF18 


C 2 + 


M 




C+C+M 


3.02xl0~ 9 





63980.7 


NIST 


TF19 


CS+ 


M 




C+S+M 


4.40XKT 10 





98600.0 


E as TF13 


TF20 


CN-i 


- M 




C+N+M 


3.32X10" 10 





74989.0 


NIST 


TF21 


SiS4 


- M 




Si+S+M 


4.40X10" 10 





98600.0 


E as TF13 


TF22 


SiN- 


h M 




Si+N+M 


3.32X10- 10 





74989.0 


E as TF20 


TF23 


S2 + 


M 




S+S+M 


7.95X10- 11 





38749.0 


NIST 


TF24 


N 2 + M 




N+N+M 


2.52X10" 7 


-1.6 


57005.4 


NIST 


NEUTRAL-NEUTRAL 


NN1 


H + 


OH 


— > 


H 2 + O 


6.86X10" 14 


2.8 


1949.5 


NIST 


NN2 


H + 


H 2 


— > 


OH+ H 2 


6.82xl0~ 12 


1.6 


9720 


NIST 


NN3 


H + 


2 


— > 


OH+ O 


6.73X10- 10 


-0.6 


8151.0 


NIST 


NN4 


H + 


CH 


— > 


C + H 2 


1.31xl0- 10 


-1.6 


80.0 


NIST 


NN5 


H + 


C 2 


— > 


C+CH 


4.67xl(r 1() 


0.5 


30450.0 


UDFA06 


NN6 


H + 


CH2 


— ► 


CH + H 2 


LOOxlO" 11 





-899.6 


NIST 


NN7 


H + 


C 2 H 


— ► 


C 2 + H 2 


5.99xl0~ n 





11800.0 


NIST 


NN8 


H + 


C2H2 


— > 


C 2 H + H 2 


lxlO" 10 





14000.0 


NIST 


NN9 


H + 


HC2O 




O+C2H2 


2.51X10- 10 








NIST 


NN10 


H + 


HC2O 




CO+CH2 


4.98X10" 11 








NIST 


NN11 


H + 


HCN 




CN+ H 2 


6.19X10" 10 





12499.0 


NIST 


NN12 


H + 


CO 




C+OH 


l.lXlO" 10 


0.5 


77700.0 


UDFA06 


NN13 


H + 


CO 




CO+H 


l.OxlO -16 








E 


NN14 


H + 


C0 2 




CO+OH 


2.51X10- 10 





13349.0 


NIST 


NN15 


H + 


cs 




CH+S 


1.2X10" 11 


0.6 


5880.0 


UDFA06 


NN16 


H + 


SiH 




Si + H 2 


2.00X10" 11 








WC98 


NN17 


H + 


SiO 




SiH + O 


LOOxlO" 16 








E as NN13 


NN18 


H + 


SiO 




Si + OH 2 


l.OOxlO" 16 








NN13 


NN19 


H + 


SiS 




SiH+ S 


l.OOxlO" 16 








NN13 


NN20 


H + 


so 




OH+ S 


5.90xl(r 1() 


-0.3 


11100 


UDFA06 


NN21 


H + 


so 




OH+ S 


5.90X10" 10 


-0.3 


11100 


UDFA06 


NN22 


H + 


NO 




OH+ N 


2.82X10- 10 





24560.0 


NIST 


NN23 


O + 


H 2 




OH+H 


3.39xl0~ 13 


2.7 


3159.3 


NIST 


NN24 


O + 


OH 




2 + H 


4.55xl0~ 12 


0.4 


-371.6 


NIST 


NN25 


O + 


H 2 




OH+ OH 


1.85X10- 11 


1 


8750 


NIST 


NN26 


O + 


CH 




OH+ C 


2.52X10" 11 





2381.0 


NIST 


NN27 


0+ CH 




CO+ H 


6.59X10" 11 








NIST 


NN28 


0+ CH 2 




CO+ H 2 


1.33X10" 11 








RH01 


NN29 


0+ CH 2 




CO+ H+H 


1.33X10 -10 








NIST 


NN30 


O + 


C 2 H 




CO+ CH 


1.69X10- 11 








NIST 


NN31 


0+ C2H2 




CO+ CH 2 


2.66X10" 10 





1980.1 


NIST 


NN32 


0+ C 2 H 2 




HC2O+ H 2 


1.19X10" 12 


2.0 


956.1 


NIST 


NN33 


0+ HC2O 




C 2 H+ 2 


1.66X10" 10 








NIST 


NN34 


0+ HC2O 




CO2H+ CH 


1.99X10" 12 








NIST 


NN35 


O + 


HCN 




OH+ CN 


1.43xl(r 12 


1.5 


3799.1 


NIST 


NN36 


O + 


SiH 




SiO+ H 


l.OOxlO -10 








UDFA06 



-50- 



Table 9 — Continued 







Chemical processes 




A 


V 




Reference a 


TV T TV T O T 

NN37 


- 


j- uu 




2 + 


C 


l.OOx 10 6 








1 — ', A TTV T -1 O 

E as NN13 


NN38 


- 


r 0U2 




C04 


o 2 


2.81xl0~ 11 





24458.0 


NIST 


NN39 


- 






C04 


■ c 


5.99xl0 -10 








NIST 


NN40 


- 


\- Ob 


* 


C04 


- s 


5.00x 10~ n 








NIST 


NN41 


- 






S04- 


c 


4.68x 10 -11 


0.5 


28940.0 


UDFA06 


NN42 


- 






C04 


- N 


1.45xl0 -10 


-0.18 





NIST 


NN43 


- 


H blO 




2 4- 


Si 


l.OOx 10 -16 








E as NN13 


TV TAT A A 

NN44 


- 


4 bllN 




Si04 


- N 


6.64x 10~ n 








NIST 


NN45 


- 


h bllN 




N04 


- Si 


5.00x 10 -11 








UDFA06 


NN46 


- 


r- bU 




2 4- 


S 


2.13x 10~ 13 


1.5 


2536.4 


NIST 


NN47 


- 


I G G 
r b2b 




S04- 


S 


1.70x 10 -11 








UDFA06 


TV T TV T A CI 

NN48 


- 


(- IN W 




2 4- 


N 


2.74x 10~ n 





21286.0 


NIST 


TV T TVT A f\ 

NN49 


- 


1VT 




N04- N 


3.01xl0 -10 


1.0 


38244.0 


NIST 


NN50 


- 


QiTvT 
4 bllN 




Si04 


- N 


6.64x 10 -11 








NIST 


NN51 


c - 






CH4 


H 


6.64xl0 -10 





11699.2 


NIST 


NN52 


c - 


h Uri 




C04 


II 


l.lOx 10~ n 








UDFA06 


NN53 


c - 






CH4- 


O 


2.25x 10 -11 


0.5 


14800.0 


UDFA06 


NN54 


c - 


1 u c \ 




CH4- OH 


1.20xl0 -12 





19776.0 


NIST 


NN55 


c - 


r Oil2 




CH4- 


CH 


2.69x 10~ 12 


0.5 





NIST 


NN56 


c - 


r U2 




C04 


O 


2.46xl0 -12 


1.5 


-613.0 


UDFA06 


NN57 


c - 






C 2 4 


- O 


5.43x10 ' 


-1.5 


57200.0 


Hanson 1973 


NN58 


c - 


r OU2 




C04 


CO 


l.lOxlO -15 








NIST 


NN59 


c - 






C04 


■ Si 


l.OOxlO -16 








1 — ', tv T TV T -1 O 

E as NN13 


TV T TVT/"* r\ 

NN60 


c - 


h bU 




C04 


- s 


3.50x 10 -11 








UDFA06 


NN61 


c - 






CS4- 





3.50x 10~ n 








UDFA06 


TV T TVT/" 1 C\ 


c - 


G 

r S2 




CS4- 


S 


7.11xl0~ n 








UDFA06 


NN63 


c - 


G~ 

h b2 




CS4- 


s 


7.11xl0~ n 








UDFA06 


NN64 


c - 






C04 


- N 


3.49x 10~ n 








NIST 


NN65 


c - 


i Mn 
h IN w 




CN4 


O 


5.57x10 


-0.3 





NIST 


TV T TVT/"* f~> 


c - 


h N 2 




CN4 


■ N 


8.69x 10 -11 





22600.0 


NIST 


TV T TVT r* r-r 

NN67 


Si - 


f H 2 




SiH4 


- II 


6.64xl0 -10 





11699.0 


E as NN51 


NN68 


Si - 


f OH 




Si04 


- H 


3.32xl0~ 10 








NIST 


NN69 


Si - 


f o 2 




Si04 


- O 


1.72X10" 10 


-0.5 


16.8 


NIST 


NN70 


Si - 


f CO 




Si04 


- C 


1.30xl0~ 9 





34516.0 


NIST 


NN71 


Si - 


f C0 2 




Si04 


- CO 


9.96xl0~ 10 





9419.1 


NIST 


NN72 


Si - 


f s 2 




SiS4- 


s 


7.00X10" 11 








E as NN62 


NN73 


Si - 


f NO 




Si04 


- N 


5. 31X10" 11 





1779.9 


NIST 


NN74 


Si - 


f O2 




Si04 


- O 


1.72X10" 10 


-0.5 


16.8 


NIST 


NN75 


S 4 


- OH 




S04- 


H 


6.59X10- 11 








NIST 


NN76 


S 4 


- O2 




S04- 


O 


1.39xl0~ 12 


0.5 


15.0 


NIST 


NN77 


S 4 


- CH 




CS4- 


II 


6.59X10" 11 








E as NN27 


NN78 


S 4 


- CO 




S04- 


C 


1.00xl0~ 16 








E as NN13 


NN79 


S 4 


- CO 




CS4- 


O 


l.OOxlO" 16 








E as NN13 


NN80 


S 4 


- cs 




S 2 + c 


1.73X10" 11 


0.5 


11500 


E as NN83 


NN81 


S 4 


- CN 




CS4- 


N 


1.45X10" 10 


-0.2 





E as NN42 


NN82 


S 4 


- SiS 




S24- Si 


4.68X10" 11 


0.5 


22 


E as NN44 
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Table 9 — Continued 







Chemical 


processes 


A 




V 


E a 


Reference a 


XTTVTO O 

JN1N83 


S + 


SO 


-> S 2 4~ O 


1.73x10" 


11 


0.5 


11500.0 




T TTlT? A flfi 

U DrAUu 


NN84 


S + 


NO 


-> S04- N 


1.75x10" 


10 





20200.0 




U DrAUu 


NN85 


N + OH 


-+ NO4- H 


4.00x 10" 


11 










INIO 1 


NN86 


N + HCN 


-> N 2 4- CH 


3.00x 10" 


10 








rL 


as ININyi 


NN87 


N + 2 


-> N04- O 


3.44x 10" 


12 


1.2 


4000 




INlo 1 


NN88 


N + CO 


-» NO4- C 


3.84x10 


_g 





35959.0 




INlo 1 


NN89 


N + CO 


-> CN4- O 


3.83x10 


_9 





35990.0 




INlo 1 


NN90 


N + C0 2 


TvT /~\ 1 f~*i /"\ 

-> N04- CO 


3.20x 10" 


13 





1710.5 




INIO 1 


NN91 


N + CN 


-* N 2 4- C 


3.00x 10" 


10 










INlo 1 


NN92 


N + CS 


-> CN4- S 


3.80x 10" 


11 


0.5 


1160.0 




T TTM? A C\£i 


NN93 


N 4 


- SiO 


-> SiN4- O 


3.84x10 


— <) 





35959.0 


i? 

rL 


as In in 00 


NN94 


N + SiO 


-> N04- Si 


l.OOx 10" 


10 








S2J 


_ _ T\TTVT1 Q 

as IN IN lo 


JN1N94 


N + SiS 


-» SiN4- S 


3.80X 10 


11 





1160.0 


li 


as INlNyZ 


NN94 


N + SO 


-> N04- S 


1.73x 10" 


11 


0.5 


750 




T TTlT? A Cld 

U DrAUu 


NN95 


N + NO 


-» N 2 4- O 


7.11x10" 


11 





786.1 




INlo 1 


NN96 


H 2 


4- OH 


-> H 2 04- H 


2.97x 10" 


12 


1.2 


2370.4 




INlo 1 


NN97 


H 2 


4- CH 


-» CH 2 4- H 


1.48x 10" 


11 


1.8 


839.4 




INlo 1 


JNN98 


H 2 


+ C 2 H 


-* C 2 H 2 4- H 


8.95X 10 


13 


2.6 


129.9 




lNTTQT 
INlo 1 


TV T TV T /"\ 

NN99 


H 2 


+ co 2 


-* CH 2 4- 2 


l.OOx 10" 


10 








IP 

Ej 


_ TvTTVT-1 Q 

as IN IN lo 


TvT TvT 1 r\r\ 

NN 100 


H 2 


4- CO 


-» CH 2 4- O 


"1 r\ v > 1 /"I — 

l.OOx 10 


33 










w 


TvTTvT 1m 

NN 101 


H 2 


+ o 2 


-> OH4- OH 


3.16x 10" 


10 





21890.0 




INlo 1 


NN 102 


H 2 


+ C 2 


-» C 2 H4- H 


l.lOx 10" 


10 





4000.0 




INlo 1 


NN103 


H 2 


4- CN 


T T /~1 TvT 1 T T 

-» HCN4- H 


5.65x 10" 


13 


2. 1 


1129.9 




KTTQT 
IN IO 1 


NN104 


OH 


4- CH 


-> H 2 04- C 


l.OOx 10" 


10 








w 
rj 


„ r1 TvTTvT-1 O 

as IN IN lo 


TvT TvT 1 r\ f\ 

NN 109 


OH 


+ C 2 H 2 - 


-» HC 2 04- H 2 


1.91x 10" 


13 





-0 




INlO 1 


TvTTvT mr 

NN 105 


OH 


+ OH 


-> H 2 4- 2 


1.65x 10" 


12 


1.1 


6013.2 




U UrAUD 


NN106 


OH 


+ OH 


-> H 2 0+ 2 


1.65x 10" 


12 


1.1 


50.0 




U DrAUu 


TvT TvT 1 r\ r ~7 

NN 107 


OH 


4- CO 


-> H4- C0 2 


3.75x 10" 


14 


1.6 


-401.9 




lNTTQT 
IN IO 1 


TvT TvT 1 AO 

NN 108 


OH 


+ CO 


-* CH4- 2 


l.OOx 10 


10 








IV 

rj 


„ rl M1NT1 Q 

as IN IN lo 


TVTTVT 1 1 f\ 

NN 110 


OH 


+ CN 


-» H+ C0 2 


3.75x 10" 


14 


1.6 


-401.9 




NIST 


NN111 


OH 


4- HCN — 


-> H 2 04-CN 


2.41x10" 


11 





5499.7 




NIST 


NN112 


o 2 


4- CH 


-» C04- OH 


8.00x10" 


11 










NIST 


NN113 


o 2 


+ CH 2 


-4 H 2 4- C0 2 


1.33x10" 


11 










RH01 


NN114 


2 + CH 2 


-» H 2 04- CO 


2.54x10" 


10 


-3.3 


1443.0 




RH01 


NN115 


o 2 


+ CH 2 


-> C0 2 4- H+H 


1.33x10" 


11 










RH01 


NN116 


o 2 


+ CH 2 


-> C04- OH4-H 


1.33x10" 


11 










RH01 


NN117 


o 2 


+ C 2 H 


-> HC 2 OH- O 


1.00x10" 


12 










NIST 


NN118 


o 2 


4- CO 


-> C0 2 4- 2 


4.20x10" 


12 





24053.0 




NIST 


NN119 


o 2 


+ N 2 


-» N04- NO 


1.00x10" 


16 








E 


as NN13 


NN120 


CH : 


2 + CH 2 - 


-4 C 2 H 2 4-H+H 


3.32x10" 


10 





5229.8 




NIST 


NN121 


CH; 


! + CH 2 - 


-> C 2 H 2 4-H 2 


2.62x10 


-9 





6009.6 




NIST 


NN122 


CO 


4- CH 


-> C 2 H4-0 


1.00x10" 


16 








E 


as NN13 


NN123 


CO 


4- CH 2 - 


-> C 2 H 2 +0 


1.00x10" 


10 










NIST 


NN124 


CO 


4- SiO 


-> C0 2 4-Si 


1.00x10" 


16 








E 


as NN13 


NN125 


CO 


4- NO 


-* C0 2 4-N 


1.00x10" 


10 








E 


as NN13 


NN126 


CN 


4- H 2 - 


-» HCN4- OH 


1.25x10" 


11 





3719.0 




NIST 
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Table 9 — Continued 



Chemical processes 



Reference 



NN127 


NO 


f NO 




N 2 + 


o 2 


2.51XKT 11 





30653.0 


UDFA06 


RADIATION ASSOCIATION 


RA1 


H + H+ 




HeH^ 


■ 4- hv 


5.26XKT 20 


-0.5 





UDFA06 


RA2 


He 4 


- H+ 




H+ 4 


- hv 


5.13xl(r 19 


1.8 





UDFA06 


RA3 


O + 







o 2 4- 


hv 


l.ooxio- 19 








Babb 1995 


RA4 


O + 


c 




CO 4 


- hv 


1.58X10" 17 


0.3 


1297.4 


Dalgarno 1990 


RA5 


O + 


c+ 




CO+ 


+ hv 


3.14X10" 18 


-0.1 


68.0 


Dalgarno 1990 


RA6 


O + 


Si 




SiO H 


- hv 


5.52X10- 18 


0.3 





Andreazza 1995 


RA7 


O + 


Si+ 




SiO+ 


+ hv 


5.50X10" 18 








Andreazza 1995 


RA8 


O + 


s 




SO 4- 


hv 


l.llxlO" 19 


0.3 


1297.9 


Andreazza 2005 


RA9 


C + 


o+ 




CO+ 


+ hv 


3.14xl0^ 18 


-0.1 


68.0 


UDFA06 


RA10 


c + 


c 




c 2 + 


hu 


4.36X10" 18 


0.3 


161.3 


Singh 2000 


RAll 


c + 


c 2 




c 3 + 


hu 


l.OOxlO^ 17 








Clayton 1999 


RA12 


c + 


s 




cs + 


hv 


4.36XKT 19 


0.2 





Andreazza 2005 


RA13 


c + 


N 




CN 4 


- hv 


7.87X10- 19 





96.0 


Singh 2000 


RA14 


Si + 


s 




CS + 


hv 


1.05X10" 19 


0.3 


66.1 


Andreazza 2007 


RA15 


s + 


S 




S 2 + 


hv 


1.38xl0^ 19 


0.3 


-78.8 


Andreazza 2005 


ION-MOLECULE 


IM1 


H 2 + - 


f He 




HeH^ 


■ + H 


1.30X10 -10 








UDFA06 


IM2 


H 2 + - 


f H 2 




H+ 4 


- H 


2.08X10" 9 








UDFA06 


IM3 


H 3 + - 


f O 




OH+ 


4- H 2 


8.40X10" 10 








UDFA06 


IM4 


H+ + Mg 




Mg+ 


4- H 2 4- H 


l.ooxio- 9 








UDFA06 


IM5 


Hj- 


f Fc 




fe+ 4 


- H 2 + H 


4.90xl0~ 9 








UDFA06 


IM6 


HcH 


+ +H 




H+ 4 


- He 


9.10X10" 10 








UDFA06 


IM7 


HeH 


+ + H 2 




H+ 4 


- He 


1.50xl0~ 9 








UDFA06 


IM8 


He+ 


+ H 2 




He 4- 


H+ 4- H 


l.lOxlO" 9 








UDFA06 


IM9 


He+ 


+ OH 




0+ 4- H 4- He 


1.10xl0~ 9 








UDFA06 


IM10 


He+ 


+ OH 




0+ H+ 4- He 


l.lOxlO -9 








UDFA06 


IM11 


He+ 


+ H 2 




OH 4 


- H+ 4- He 


2.04X10" 10 








UDFA06 


IM12 


He+ 


+ o 2 




0+ 4 


- O 4- He 


l.ioxicr 9 








UDFA06 


IM13 


He+ 


+ CO 




C+ 4 


O + He 


1.40X10" 9 


-0.5 





UDFA06 


IM14 


He+ 


+ co 2 




CO+ 


4- O 4- He 


8.70X10- 10 








UDFA06 


IM15 


He+ 


+ c 2 




C+ 4 


- C 4- He 


1.60X10" 9 








UDFA06 


IM16 


He+ 


+ CH 2 




C+ 4 


- H 2 4- He 


7.50X10- 10 








UDFA06 


IM17 


He+ 


+ C 2 H 2 




C+ 4 


- H 2 4- He 


1.61X10" 9 








UDFA06 


IM18 


He+ 


+ SiO 




0+ 4 


- Si+ He 


8.60X10" 10 








UDFA06 


IM19 


He+ 


+ SiO 




Si+ 4 


- O 4- He 


8.60X10" 10 








UDFA06 


IM20 


C+ - 


f O 2 




CO+ 


4- O 


3.80X10- 10 








UDFA06 


IM21 


C+ - 


f O 2 




0+ 4 


- CO 


6.20X10- 10 








UDFA06 


IM22 


C+ - 


f CO 2 




CO+ 


4- CO 


l.ioxicr 9 








UDFA06 


IM23 


c+ - 


f SiO 




Si+ 4 


- CO 


5.40X10- 10 








UDFA06 


IM24 


c+ - 


f SO 




S+ + 


CO 


2.60X10- 10 








UDFA06 


IM25 


c+ - 


f SO 




CO+ 


4- S 


2.60X10" 10 








UDFA06 


IM26 


Si+ 


+• OH 




SiO+ +H 


6.30X10" 10 








UDFA06 


IM27 


s+ - 


- OH 




SO+ 


4- H 


6.10X10" 10 








UDFA06 


IM28 


s+ - 


- o 2 




SO+ 


4- O 


1.50X10" 11 








UDFA06 
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Table 9 — Continued 







Chemical 


processes 


A 


V 


E a 


Reference a 


IM29 


(V 


4- O 


-> CCT 


- + c 


3.ioxicr 10 








UDFA06 


IM30 


SiO 


+ + O 


-> Si+ 


+0 2 


2.ooxio- 10 








UDFA06 


IM31 


SiO 


+ 4- C 


-> Si+ 


+ CO 


l.OOxlO" 9 








UDFA06 


IM32 


SiO 


+ 4- CO - 


- Si+ 


4- CO2 


7.90xl(r 10 








UDFA06 


IM33 


SiO 


+ 4- S 


-> Si+ 


+ SO2 


l.ooxicr 9 








UDFA06 


IM34 


SiO 


+ 4- N 


-> Si+ 


4- NO 


2.ioxicr 10 








UDFA06 


CHARGE EXCHANGE 


CE1 


H 4 


-H+ — H+ 


+ H 2 


6.40X1CT 10 








UDFA06 


CE2 


O 4 


-c 2 + 


-> co- 


- +C 


3.ioxicr 10 








UDFA06 


CE3 


O 4 


- CO+ 


-> o+ 


4- CO 


1.40X10 -10 


0.5 





UDFA06 


CE4 


O 4 


■c 2 + 


-+ CO+ 4- C 


3.10X10- 10 








UDFA06 


CE5 


C 4 


CO+ 


- c+ 


4- CO 


l.lOxlO -10 








UDFA06 


CE6 


C 4 


c 2 + 


-> c+ 


+ c 2 


l.lOxlO -10 








UDFA06 


CE7 


C 2 


+ o+ 


- c+ 


+ 


4.80X10" 10 








UDFA06 


CE8 


CO 


4-0+ 


-> CO+ + o 


4.90X10" 12 


0.5 


4580.0 


UDFA06 


CE9 


Si 4 


- C+ 


-» Si+ 


+ c 


2.10X10" 9 








UDFA06 


CE10 


Si 4 


- S+ 


-» Si+ 


4- S 


i.eoxio -9 








NIST 


CE11 


SiO 


4- H+ 


-> SiO+ + H 


3.30X10" 9 








UDFA06 


CE12 


S + 


C+ 


-* S+ - 


f C 


1.50X10" 9 








NIST 


CE13 


s + 


CO+ 


-* S+ - 


f CO 


l.lOxlO -9 








NIST 


CE14 


Fc + 0+ 


-> Fe+ 


+ 


2.90X10" 9 








UDFA06 


CE15 


Fc - 


hC+ 


-> Fe+ 


+ c 


2.60X10" 9 








UDFA06 


CE16 


Fe - 


4 Si+ 


-> Fe+ 


4- Si 


1.90xl0~ 9 








UDFA06 


CE17 


Fe - 


r s+ 


-> Fe+ 


+ S 


1.80X10- 10 








NIST 


CE18 


Fe - 


h SO+ 


-> Fe+ 


4- SO 


1.60X10" 9 








UDFA06 


CE19 


Mg 


+ H+ 


-» Mg H 


- + H 


l.lOxlO" 9 








UDFA06 


CE20 


Mg 


4-0+ 


-> Mg H 


- + O 


l.lOxlO" 9 








UDFA06 


CE21 


Mg 


+ C+ 


-> 


■ + C 


l.lOxlO -9 








UDFA06 


CE22 


Mg 


4-Si+ 


-+ Mg H 


- + Si 


2.90X10" 9 








UDFA06 


CE23 


Mg 


4-SiO+ 


-» Mg H 


- 4- SiO 


l.OOxlO -9 








UDFA06 


CE24 


Al - 


1- o+ 


-> A1+ 


+ O 


2.00X10" 9 








UDFA06 


CE25 


Al - 


1- c+ 


-* A1+ 


+ c 


2.70X10" 9 








UDFA06 


CE26 


Al - 


(- CO+ 


-* A1+ 


4- CO 


2.20xl0~ 9 








UDFA06 


CE27 


Al - 


F Si+ 


-> A1+ 


+ Si 


2.70X10" 9 








UDFA06 


CE28 


Al - 


(- SiO+ 


-> A1+ 4- SiO 


2.20X10" 9 








UDFA06 



ELECTRONIC RECOMBINATION 



ER1 


H+ + 


e — 


— > H 




3.50X10" 12 


-0.7 





UDFA06 


ER2 


H+ + 


e — 


H 4 


H 


1.60X10" 8 


1.4 





UDFA06 


ER3 


H+ + 


e — 


— » H 2 


+ H 


2.34xl0~ 8 


-0.5 





UDFA06 


ER4 


H+ + 


e — 


H 2 


+ H 


2.34XKT 8 


-0.5 





UDFA06 


ER5 


C+ + 


e — 


— > C 




4.67X10" 12 


-0.6 





UDFA06 


ER6 


ct + 


e — 


-> C 4 


C 


3.00X10" 7 


-0.5 





UDFA06 


ER7 


CO+ 


+ e- 


-> C 4 





2.00xl0~ 7 


-0.5 





UDFA06 


ER8 


0+ + 


e — 


— > O 




3.24X10" 12 


-0.7 





UDFA06 


ER9 


Si+ 4 


e — 


— > Si 




4.90X10" 12 


-0.6 





UDFA06 


ER10 


SiO+ 


+ e- 


-» Si 4 


- 


2.00X10" 7 


-0.5 





UDFA06 
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Table 9 — Continued 







Chemical processes 




A 


V 


E a 


Reference a 


ERll 


S+ + e" 


— > S 






3.90xl(r 12 


-0.6 





UDFA06 


ER12 


SO+ + c 


- — ♦ s + 







2.00xl0~ 7 


-0.5 





UDFA06 


ER13 


Fe+ + e" 


— > Fe 






3.70X10- 12 


-0.6 





UDFA06 


ER14 


Mg+ + c 


- — ► Mg 






2.80xl0~ 12 


-0.9 





UDFA06 


ER15 


A1+ + e- 


— > Al 






3.24X10" 12 


-0.7 





UDFA06 


COMPTON ELECTRON DESTRUCTION 


CED1 


O 


— ► 0+ 


f e" 




1.57X10" 5 





3464.1 


See text 


CED2 


C 


— » C+ - 


f e" 




1.99X10" 5 





3464.1 




CED3 


Si 


— > Si+ 


+ e" 




1.57X10 -5 





3464.1 




CED4 


s 


— ► s+ - 


h e - 




1.57X10" 5 





3464.1 




CED5 


Fe 


— > Fc+ 


+ e" 




1.57xl0~ 5 





3464.1 




CED6 


Mg 


— > Mg+ 


+ e- 




1.57xl0~ 5 





3464.1 




CED7 


Al 


— ^ A1+ 


+ e- 




1.57xl0" 5 





3464.1 




CED8 


H_> 


— > H + 


H+ + 


c 


8.86xl0~ 7 





3464.1 




CED9 


H 2 


— H+- 


f e- 




1.93xl0~ 5 





3464.1 




CED10 


H 2 


— > H + 


H 




9.43xl0~ 6 





3464.1 




CED11 


OH 


— > H + 


0+ + 


c~ 


9.46X10" 7 





3464.1 




CED12 


OH 


— > O + 


H- + 


c~ 


2.94xl0~ 6 





3464.1 




CED13 


OH 


— > O + 


H 




5.81xl0~ 6 





3464.1 




CED14 


H 2 


— > OH 


■f H+ ■ 


I- e- 


2.94X10" 6 





3464.1 




CED15 


H 2 


— > OH 


■f H 




5.81xl0~ 6 





3464.1 




CED16 


CH 


— > H + 


C+ + 


c~ 


9.46xl0~ 7 





3464.1 




CED17 


CH 


— > C + 


H+ + 


c~ 


2.94X10" 6 





3464.1 




CED18 


CH 


— » H + 


C 




5.81X10" 6 





3464.1 




CED19 


CH 2 


— ► CH 


f H+ - 


f e- 


2.94xl0~ 6 





3464.1 




CED20 


CH 2 


— > CH 


f H 




5.81X10" 6 





3464.1 




CED21 


C 2 H 2 


— ► C 2 H + H+ 


+ e- 


2.94X10" 6 





3464.1 




CED22 


C2H2 


— > C 2 H 


+ H 




5.81xl0~ 6 





3464.1 




CED23 


HCO 


— > CO 


-rH+ ■ 


f e- 


2.94X10" 6 





3464.1 




CED24 


HCO 


— > H + 


CO+ 


f e- 


8.86XIO- 7 





3464.1 




CED25 


HCO 


— > H + 


CO 




5.81xl0~ 6 





3464.1 




CED26 


O2 


— > O + 


0+ + 


c~ 


9.46xl0~ 7 





3464.1 




CED27 


o_. 


— > O + 







5.81X10" 6 





3464.1 




CED28 


CO 


— > C + 


0+ + e~ 


9.46xl0~ 7 





3464.1 




CED29 


CO 


— > O + 


C+ + 


c 


2.94xl0~ 6 





3464.1 




CED30 


CO 


— > CO+ + e- 




2.14xl0~ 5 





3464.1 




CED31 


CO 


— > C + 







5.81xl(r 7 





3464.1 




CED32 


co 2 


— ► CO + 0+ 


+■ e- 


9.46xl0~ 7 





3464.1 




CED33 


co 2 


— > + 


CO+ 


■f e" 


2.94X10" 6 





3464.1 




CED34 


co 2 


— ► CO + 




5.81X10" 7 





3464.1 




CED35 


co 2 


— ► CO + 0+ 


f e" 


9.46xl0~ 7 





3464.1 




CED36 


c 2 


— > c + 


c+ + 




9.46X10" 7 





3464.1 




CED37 


c 2 


- c+- 


f e- 




2.14X10" 5 





3464.1 




CED38 


c 2 


— > c + 


C 




5.81X10" 6 





3464.1 




CED39 


SiO 


— > Si + 


0+ + 


c~ 


1.07xl0~ 7 





3464.1 




CED40 


SiO 


— -> + 


Si- + 




3.33X10" 6 





3464.1 
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Chemical processes 




A 


;/ 


E a Reference a 


CED41 


SiO 


> SiO+ + e" 




2.42xl0~ 5 





3464.1 


CED42 


SiO 


► O + 


Si 




6.58X10" 6 





3464.1 


CED43 


SiS 


— ► Si 4 


S+ 4 


- e 


2.94X10 -6 





3464.1 


CED44 


SiS 


— > S + 


Si+ 4 


- e — 


2.94xl0" 6 





3464.1 


CED45 


SiS 


— > Si 4 


S 




5.81X10" 6 





3464.1 


CED46 


SO 


— > + 


S+ 4 


e — 


2.94x10-° 





3464.1 


CED47 


so 


— > s + 


0+ 4 


e — 


9.46X1CT 7 





3464.1 


CED48 


so 


— ► so 4 


' 4- e" 




2.14X10" 5 





3464.1 


CED49 


so 


— > + 


S 




5.81x10-° 





3464.1 


CED50 


so 


— > + 


S+ 4 


e — 


2.94xl0" 6 





3464.1 


CED51 


s 2 


— > s + 


S+ + 




9.46xl0" 7 





3464.1 


CED52 


S 2 


— ► s + 


S 




5.81xl0~ 6 





3464.1 


CED53 


N 2 


— > N + 


N 




5.44X10" 6 





3464.1 


UV PHOTODISSOCIATION 


PHOT1 


H 2 


— ► H+ H 




3.80X10" 5 





3464.1 See text 


PHOT2 


2 


— ► O + 


O 




3.80X10" 5 





3464.1 


PHOT3 


OH 


— > O + 


H 




3.80X10- 5 





3464.1 


PHOT4 


H 2 


— ► OH 


4- H 




3.80X10" 5 





3464.1 


PHOT5 


CO 


— ► C + 


O 




3.80X10" 5 





3464.1 


PHOT6 


co 2 


— ► CO 


4- O 




3.80X10" 5 





3464.1 


PHOT7 


c 2 


— > c + 


C 




3.80X10- 5 





3464.1 


PHOT8 


co 2 


— > CO 


+ O 




3.80xl0~ 5 





3464.1 


PHOT9 


SiO 


— ► Si 4 


O 




3.80xl0~ 5 





3464.1 


PHOT10 


so 2 


— > s + 


O 




3.80X10- 5 





3464.1 


PHOT11 


N 2 


— > N + 


N 




3.80X10" 5 





3464.1 



a NIST is the NIST chemical kinetics database. UDFA06 is by Woodall et al. (2007). 
Other references are listed in the bibliography. 'E' means 'estimated'. 



